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ABSTRACT 

We present new nonlinear, time-dependent convective hydrodynamical models of RR Lyrae stars computed 
assuming a constant helium-to-metal enrichment ratio and a broad range in metal abundances (Z = 0.0001- 
0.02). The stellar masses and luminosities adopted to construct the pulsation models were fixed according to 
detailed central He burning Horizontal Branch evolutionary models. The pulsation models cover a broad range 
in stellar luminosity and effective temperatures and the modal stability is investigated for both fundamental 
and first overtones. We predict the topology of the instability strip as a function of the metal content and new 
analytical relations for the edges of the instability strip in the observational plane. Moreover, a new analytical 
relation to constrain the pulsation mass of double pulsators as a function of the period ratio and the metal content 
is provided. We derive new Period-Radius-Metallicity relations for fundamental and first-overtone pulsators. 
They agree quite well with similar empirical and theoretical relations in the literature. From the predicted 
bolometric light curves, transformed into optical (UBVRI) and near-infrared ( JHK ) bands, we compute the 
intensity-averaged mean magnitudes along the entire pulsation cycle and, in turn, new and homogenous metal- 
dependent ( RIJHK ) Period-Luminosity relations. Moreover, we compute new dual and triple band optical, 
optical-NIR and NIR Period-Wesenheit-Metallicity relations. Interestingly, we find that the optical Period- 
W (V,B - V) is independent of the metal content and that the accuracy of individual distances is a balance 
between the adopted diagnostics and the precision of photometric and spectroscopic datasets. 

Subject headings: stars: evolution — stars: horizontal-branch — stars: oscillations — stars: variables: RR 
Lyrae 


1. INTRODUCTION 


RR Lyrae stars (RRLs) are low-mass helium burning 
stars playing a crucial role both as standard candles and 
tracers of old (t>10 Gyr) stellar populations dBono et all 


different Galactic (see e.g. Vivas & Zinn 

2006; IZinn et alJ 

2014| Drake et al. 2013 Pietrukowicz et al. 

2014, and refer- 

ences therein) and extragalactic (see e.g. Moretti et al. 2009; 

Soszvnski et al. :2009; 

Fiorentino et al. 2010 

Soszvnski et al. 

2010; Fiorentino et al 

2012; Cusano et al. 2013, and refer- 


ences therein) envir onments, including a significant fraction 
of globular clusters (Coppola et al.||201 1; |Di Criscienzo et all 


120111 iKunder et all 2013b iKuehnetalJ 2013h . One of the 
key advantage in using RRLs is that they can be eas¬ 
ily identified thanks to the shape of their light curves, 
the luminosity amplitudes and pulsation periods. They 
are also relatively bright, typically 3.0/3.5 mag brighter 
than the Main Se quence Turn Off (MSTO) stars (see 
e.g. iPiersimoni et aT]l2002! iRineni et~aTll2007b ICoppola et all 

12011b 120131 iBraga et all 12014 ). The RRLs have been pop¬ 
ular primary distance indicators thanks to the relation be¬ 
tween the absolute visual m agnit ude and the iron ^bun- 
dance [Fe/H] f Caputo et al J 200(1 ICacciari & Clementinil 
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120031 iBono et ahl 1201 lb IMarconil l2012h . The intrinsic er¬ 
rors and systematics affecting distances based on this rela- 
tion have been widely discussed in the literatu re (Cassi si et al 


_1998; [Caputo et al. 2000; Di Criscienzo, Marcon i & Canuto 
2004 : ICassisi et all 12008b IMarconil 12009b 1201 2l) However. 


RRLs have been empirically recognised to obey to a 
Period -Luminosity (PL) rela t ion in the near-infrared (NIR) 
bands (Longm ore et al.lll986L 1 1990h . The physical bases, the 
key advantages in using NIR PL relations to estimate indi¬ 
vidual distances together with their metallicity dependence 
have been extensively discussed in the literature both from the 


2001, 2002, 2001 

Catelan et al. 2004[ 

Dall’Ora et al. 2004 

Sollima et al. 2006> 

Marconi 2009, 2012 

; Coppola et al. 2011 


Bono et all 1201 lb IBraga et aP 120141) . Here, we only men¬ 

tion that the NIR PL relations are marginally affected by 
uncertainties in reddening correction and by evolutionary 
effects (IBono et all l2003ii . In spite of these unquestion¬ 
able advantages the NIR PL relations might also be prone 
to systematic errors. a) The reddening correction be¬ 
comes a thorny problem if the targets are affected by dif¬ 
ferential reddening. This is the typical problem in deal¬ 
ing with RRL dista nces in the low-redd ening regions of 
the Ga lactic Bul ge (Matsunag a et all [201 3h and in the inner 
Bulge (Soszvn ski et all 2014) . b) Even if the width in mag¬ 
nitude of the RRL instability strip at fixed period is almost 
halved when moving from the I-band to the K-band, the PL re¬ 
lations are intrinsically a statistical diagnostics to estimate the 
distances, since the width in temperature is neglected. This 
means that the intrinsic dispersion of the PL relations, even in 
the NIR bands, is still affected by the width in temperature of 
the instability strip. 

To overcome the above problems it has been empirically 
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suggested to use the optical and the optical-NIR Period- 
Wesenheit (PW) relation s dDi Criscienzo. Marconi & Caputo I 
2004 [Braga et al] 120141 Coppola et al. 2015, in prepara¬ 
tion). The recent literature concerning the use of the red¬ 
dening free Wesenheit magnitudes is quite extensive, but it 
is mainly focusse d on classical Cepheids (iRjpepijTal j |20j2 
iRiess et afl 1201 2t iFiorentino et al.l 1201 3h ifnno et alJ 120131) . 
The two main advantages in using the PW relations are: a) 
they are reddening free by construction; b) they mimic a 
Period-Luminosity-Color (PLC) relation. This means that in¬ 
dividual RRL distances can be estimated with high accuracy, 
since they account for the position of the object inside the 
instability strip. Therefore, the PW relations have several ad¬ 
vantages when compared with classical distance diagnostics. 
However, their use to estimate distances of field and clus¬ 
ter RRLs has_been quite limited. They have been adopted 
by |Di Criscienzo. Marconi & Caputo I (120041) : iBraea et al] 
(120141) to estimate the distance of a number of Galactic G Cs 
and by ISoszvnski et alJ (120141) : iPietrukowicz et all (120141) to 
derive the distances of RRLs in the Galactic Bulge. The 
inferred individual and average distances, based on the two 
quoted methods (PL and PW relations), can be used to de¬ 
rive the 3D structure of the investigated stellar systems, as 
well as to trace r adial trends across the Halo and tidal stel¬ 
lar streams (e.g. [S oszyfekietaLl |2010|; |Cusano_elal] |2013|; 
Moretti et al.ll2014t IPietrukowicz et al.ll20l4ISoszvhski et all 

20141) . 


The lack of detailed investigations concerning pros and 
cons of RRL PW relations also applies to theory. Indeed, we 
still lack detailed constraints on their metallicity dependence 
and on their intrinsic dispersion when moving from optical, to 
optical-NIR and to NIR PW relations. 

The above limitations became even more compelling, dur¬ 
ing the last few years, thanks to several ongoing large- 
scale, long-term sky variability surveys. The OGLE IV col¬ 
laboration already released VI -band photometry for more 
than 38,000 RRLs in the Galactic Bulge (IPietrukowicz et afl 
12014 ISoszvnski et alJ 120141) . They plan to release simi¬ 
lar data for RRLs in the Magellanic Clouds (MCs). The 
ASAS collaboration is still collecting VI -band photometry 
for the entire southern sky ( iPoimariskil 120141) . Multi-band 
photometry for a huge number (~23000) of southern field 
RRLs has also been relea sed by CATALINA (Dra ke et alJ 
l2013l:lTbrrealba et al.ll2014l) . Similar findings have been pro¬ 
vided in t he NIR bands by V W for RRLs in the Galac¬ 
tic bulge dMinniti et al.l 120141) and by VMC (Moretti et al. 
2015, in preparation) for RRL in the MCs. New optical cat¬ 
alogs have also been released by large optical surveys such 
as SDSS and Pan-STARRSl (Abbas~et alj2014 l) and UV sur¬ 
veys such as GALEX dGezari et al. 11201 3t iKinman & Brownl 
120141) . New detections of field RRLs have also been provided, 
as ancillary results, by photometric surveys interested in the 
identificatio n either of moving objec ts (LINEAR, iSesar et alJ 
1201 Ul20l3l) or of near earth objects dMiceli et al.ll2008l) or of 
transient phenomena (ROTSE, Kinemuchi et al. 2006; PTF, 
Sesar et al. 2014) or of the identification of stellar stre ams in 
the Galactic Halo (IVivas & Zinnll2006l : IZinn et al.ll2014l) 

The above evidence indicates that a comprehensive theo¬ 
retical investigation addressing PL and PW relation in optical 
and NIR bands is required. Our group, during the last twenty 
years, constructed a detailed evolutionary and pulsation 


Bono et al. 1995c [Bono et aLj 1997 

Marconi et al. 

2003; 

Di Criscienzo, Marconi & Caputo 

2004; Marconi et al. 


1201 ll) . We have computed several grids of Horizontal 
Branch (HB) models and pulsation models accounting for 
a wide range of chemical compositions (iron, helium), 
stellar masses and luminosity levels. The above theoretical 


framework was ad 
properties of RRLs 
globular clusters (C 
dwarf galaxies (e. 

opted to compare predicted 
in a wide range of stellar 
iCs) (o) Cen, Marconi et al. 

and observed 
environments: 
20ll[), nearbv 

g. Carina and Hercules; 

Coppola et al. 

2013; Stetson et al. 2014 Musellaetal. 2012), the Galactic 

bulge (Bonoetal. 

19971 Groenewegen et al.l 2008) and the 

Galactic halo (IFiorentino et al. 20141). 


However, the quoted theoretical framework was built using 
a broad range of evolutionary prescriptions concerning stel¬ 
lar masses, lumin osity levels and the i r dependence on chem¬ 
ical compos ition (ICassisi et al. 1998, 2004; Pietr inferni et al.l 
20041 120061) . Therefore, we decided to provide a new spin on 
the evolutionary and pulsation properties of RRLs covering 
simultaneously a broad range in chemical compositions, stel¬ 
lar masses and luminosity levels. The current approach, when 
compared with the quoted pulsation investigations, has sev¬ 
eral differences. We take account of seven (^-enhanced chem¬ 
ical compositions ranging from very metal-poor (Z = 0.0001) 
to the standard solar value (Z = 0.02). The adopted stellar 
masses and luminosity levels, at fixed chemical composition, 
are only based on evolutionary prescriptions, rather than rely¬ 
ing on a grid of masses and luminosities encompassing evo¬ 
lutionary values. 

The structure of the paper is the following. In §2 we present 
the new theoretical framework and discuss the input physics 
adopted to construct both evolutionary and pulsational mod¬ 
els. §3 deals with the new pulsation relations for fundamental 
(FU) and first overtone (FO) pulsators. The topology of the in¬ 
stability strip as a function of the chemical composition is dis¬ 
cussed in §4. In this section we also address the role of the so- 
called ”OR region” to take account of the pulsation properties 
of mixed-mode pulsators. The predicted period-radius (PR) 
relations and the comparison with similar predicted and em¬ 
pirical relations is given in §5. In §6 we discuss the predicted 
light curves and their transformation into optical and NIR ob¬ 
servational planes. The new metallicity dependent PL and PW 
relations are discussed in §7. The summary of the results of 
the current investigation are outlined in §8 together with the 
conclusions and the future developments of this project. 

2. EVOLUTIONARY AND PULSATION FRAMEWORK 
2.1. HB evolutionary models 

The theoretical framework to compute HB evolu- 
tio nary models has alr eady been discussed in detail 
by iPietrinferni et all (120061) . The interested reader is referred 
to the above paper for a thorough discussion concerning the 
input physics adopted for central helium burning evolution¬ 
ary phases. The entire set of HB models are available in the 
BaSTI database 6 . The HB evolutionary models were com¬ 
puted, for each assumed chemical composition, using a fixed 
core mass and envelope chemical profile. They were com¬ 
puted evolving a progenitor form the pre-main sequence to 
the tip of the RGB with an age of 13 Gyr. The RGB pro¬ 
genitor has typically a mass of the order of to 0.8 M e in 
the very metal-poor regime increasing up to ~ 1.0 M e in the 
more metal-rich regime. The mass distribution of HB models 
ranges from the mass of the progenitors (coolest HB models) 

6 http://www.oa-teramo.inaf.it/BASTI/ 
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down to a total mass of the order of 0.5 M e (hottest HB mod¬ 
els). The trend in He core and envelope mass as a function 
of metal and helium abundances will be discussed in a forth¬ 
coming paper (Castellani et al. 2015, in preparation). 

The evolutionary phases off the Zero-Age-Horizontal 
Branch (ZAHB) have been extended either to the onset of 
thermal pulses for more massive models or until the luminos¬ 
ity of the model, along the white dwarf cooling sequence be¬ 
comes, for less massive structures, fainter than log (L/L G )— 
2.5. The adopt ed g-enhanced ch emical mixture is given in 
Table 1 of lPietrinferni et all (120061). The g-elements were en¬ 
hanced with respect to the iGrevesse et al.1 (119931) solar metal 
distribution by variable factors. We mainly followed elemen¬ 
tal _abundances for field old, low-mass stars by iRvan et all 
(1199 ll) . The overall enhancement—[g/Fe]—is equal to 0.4. 

To constrain the metallicity dependence of RRL pulsation 
properties we adopted seven different chemical compositions, 
namely Z= 0.0001, 0.0003, 0.0006, 0.001, 0.004, 0.008 and 
0.0198. We also assumed, according to recent CMB ex- 
periment s (lAde etliD 120141) . a primordial He-abundance of 
0.245 (Cassisi et al.ll2003t) . together with a helium-to-metals 
enrichment ratio of AY/AZ= 1.4. The adopted AY/AZ value 
allows us to match the calibra ted initial He abundanc e of the 
Sun at solar metal abundance (tSerenelli & Basull2010l) . 

Figure [I] shows the Hertzsprung-Russell diagram for three 
sets of HB models. From left to right the black solid line 
shows the location of the ZAHB, while the dashed black line 
the central helium exhaustion. The red solid lines display 
three evolutionary models of HB structures populating the 
RRL instability strip. They range from 0.75, 0.76, 0.77 M Q 
for the most metal-poor chemical composition (right panel), 
to 0.65,0.66,0.67 M 0 forZ=O.OOl (middle panel) and to 0.56, 
0.57, 0.58 for Z=0.008 (left panel). 

The above evolutionary prescriptions show three features 
relevant for RRL properties. 

a) - An increase in the metal content from Z—0.000 1 to 
Z=0.008 causes on average, in the middle of the instability 
strip (log 7’ e yy=3.83-3.85), a decrease of -0.2 dex in the mean 
luminosity level log (L/L Q ) — 1.8 and 1.6), but also a decrease 
in stellar mass (from 0.76 to 0.57 M/M a ). This well known 
evolutionary evidence (ICastellani et al.lM 99 1) and the coeffi¬ 
cients of both stellar luminosity and stellar mass in the pulsa¬ 
tion relations (see § 3 and equations 1) explain the observed 
decrease in period when moving from metalpoor to metalrich 
stellar structures dFiorentino et al.ll2014b . 

b) - An increase in the metal content from Z—0.000 1 to 
Z—0.008 causes a steady increase in the luminosity width of 
the HB. The increase is almost a factor of two in the visual 
band (~0.5 vs -1 m ag). This evidence was originally sug¬ 
gested by iSandagel (119931) and later supported by theoreti¬ 
cal models TBono, Marconi. & S te11 i n g werfl11999t) . The con¬ 
sequence of this evolutionary feature is a steady increase in 
the intrinsic spread in luminosity during central helium burn¬ 
ing phases when moving from metal-poor to the metal-rich 
structures. 

c) - An increase in the metal content from Z=0.0001 to 
Z—0.008 causes a steady increase in the width in temperature 
of the hook performed by evolutionary models in the early 
off-ZAHB evolution. The extent in temperature is ma inly 
driven by the efficiency of the H-shell burning dBono et al.1 
1 1 997t ICassisi et afllT9 98 ). The above evolutionary feature af¬ 
fects the lifetime that the different stellar structures spend in¬ 
side the instability strip, and in t urn, it opens the path to the 
so-called hysteresis mechanism (Ivan Albada & Baker! 11973b 


iBono et ak~lll995t IBono et aklll995L 119971k To constrain on 
a more quantitative basis the above effect, we estimated the 
typical central He burning time of a metal-poor (Z = 0.0001, 
0.76 M/M q ), a metal-intermediate (Z = 0.001, 0.66 M/M 0 ) 
and a more metal-rich (Z = 0.008, 0.57 M/M Q ) stellar struc¬ 
ture located in the instability strip (see column 4 in Table ??). 
We found that the He-burning lifetime steadily increases as 
a function of the metal content from -67 Myr, to -77 Myr 
and to -86 Myr. This means an increase on average of the 
order of 25% when moving from metal-poor to metal-rich 
RRLs. To constrain the difference in in RRL production rate 
when moving from a metal-poor to a metal-rich stellar pop¬ 
ulations the He-burning lifetimes need to be normalised to a 
solid evolutionary clock. We estimate the central H burning 
phase as the difference between the main sequence turn-off 
(MSTO) and a point along the main sequence that is 0.25 dex 
fainter. To estimate the H-burning lifetimes we selected the 
stellar mass at the MSTO of a 13 Gyr cluster isochrone for 
the three quoted chemical compositions. The stellar masses 
and the H-burning lifetimes are also listed in Table ??. The 
ratio between He and H evolutionary lifetimes ranges from 
0.032 for the metal-poor, to 0.023 for the metal-intermediate 
to 0.014 for the metal-rich structures. This means that on av¬ 
erage the number of RRL per MS star in a metal-poor stellar 
population is at least a factor of two larger than in a metal-rich 
stellar population. This evidence is further suggesting that the 
number of RRLs and their distribution across the instability 
strip is affected by different parameters: the topology of the 
instability strip, the excursion in temperature of HB evolution¬ 
ary models and the evolutionary lifetimes (IBono et al .11 19971 
iMarconi et al.ll201lt iFiorentino et al.ll2014l). 

2.2. Stellar pulsation models 

The pulsation properties of RRLs have been extensively 
inv estig a ted bv several authors, since the pioneering papers 
by |Co2 (179631) : Casto r {Up); van Albada & Baker (179771) 
based on linear non adiabatic models. The advent of non¬ 
linear hydrocodes provided the opportunity to investigate the 
limiting cycle behaviour of radial variables (IChristvl 1 1967k 
ICox!ll974l) . However, a detailed comparison between theory 
and observations was only possible with hydrocodes taking 
account of the coupling between convection motions and ra¬ 
dial displacements. This approach provides solid predictions 
not only on the red edges of the instability_strip but also on 
pulsation amplitudes and m o dal stability (IStelling werfll 19821: 
Bono & Stellingwerfj 119941 iFeuchtingerl 1 1 999t IBono et al.1 
2000 ). 

Recent developments in the field of nonlinear mod- 
ell ing of RRLs were provided by ISzabo et aD (120041) and 
by Sm olec et al.l (120131) . They also take account of evolution¬ 
ary and pulsational properties of RRLs using an amplitude 
equation formalism. 

The key advantage of our approach is that we simul¬ 
taneously solve the hydrodynamic al conservations equa- 


tions together with a nonlocal (mixing-length like), time- 
dependent treatment of convection transport (Stellingweii 

1982 

; Bono & Stellingwerl 1994j Bono et al. 2000; Marconi 

2009 

of RI 
in st< 
tion 

). We performed extensive and detailed computations 
^L nonlinear convective models covering a broad range 

11 ar masses, luminosity levels and chemical composi- 
Bono et al. 1997, 1998; Marconi et al. 2003, 2011). To 


constrain the impact that the adopted treatment of the convec¬ 
tive transport has on pulsation observables we have also com¬ 
puted different sets of models changing the efficiency of con- 




















































































4 


vection (iDi Criscienzo. Marconi & Caputo I 120041: iMarconil 
120091) . 

The RRL pulsation models presented in this paper 
were constructed using the hydrodynamical code developed 
bv IS tellin g werll (1 1 9821) and u pdated by |Bono_^Stellingwerl 
(1994): iBono et al.l (119981) : [Bono, Marcon i. & Stellin gwerf 
(1999, see also Smolec & Moskalik 2010 for a simi¬ 
lar approach). The physical and numerical assumptions 
adopted to compute these models are the sam e discussed 
inlBono et al.1 (ll998l):lBono. Marconi. & Stelhngwerll (119991) : 
IMarconi et al.1 (120031 1201 ilk In parti cular, we adopted 
the OP AL radiative opacities released by Iglesi as~& Rogers 1 
(1996, http://www-phys.llnl. gov/Research/OPAL/opal.html ) 
and the molecular opacities bv lAlexander & Ferguson! (119941) . 

To compute the new models we adopted seven different 
metallicities ranging from the very metal-poor regime (Z = 
0.0001) to the canonical solar abundance (Z = 0.02). The 
helium abundance, for each chemical composition, was fixed 
according to the helium-to-metals enrichment ratio adopted 
by Pietrinferni et al.1 (120061) . namely AY/AZ ~ 1.4 for a pri¬ 
mordial helium content of ~ 0.245. The new RRL models 
when compared with similar models computed by our group 
present several differences, a) The stellar mass, at fixed chem¬ 
ical composition, is the mass of the ZAHB ( Mzahb ) predicted 
by HB evolutionary models at the center of the instability strip 
(log T e ~ 3.85). The current predictions indicate that the cen¬ 
ter of the instability strip is located at log T e ~ 3.82. We 
adopted the former value, since this is the canonical value 
and makes more solid the comparison with predictions avail¬ 
able in the literature. Note that the difference on Mzahb is 
minimal when moving from log T e ~ 3.85 to log T e ~ 3.82. 
b) The faintest luminosity level, at fixed chemical composi¬ 
tion, is the predicted ZAHB luminosity level, c) To take ac¬ 
count of the intrinsic width in magnitude ( Sandage i 1990 ) of 
the RRL region we also adopt, at fixed chemical composi¬ 
tion and ZAHB mass value, a luminosity level that is 0.1 dex 
brighter than the ZAHB luminosity level, d) Theory and ob¬ 
servations indicate that stellar structures located close to the 
blue edge of the instability strip cross the instability strip dur¬ 
ing their off ZAHB evolution. These stellar structures typi¬ 
cally evolve from hotter to cooler effective temperatures and 
cross the instabilit y strip at luminosity levels higher than typ¬ 
ical RRL (a> Cen. [Marconi et akll201 ll) . To take account for 
these evolved RRLs we also adopted, at fixed chemical com¬ 
position, a second value of the stellar mass— M evo —that is 
10% smaller than Mzahb- The decrease in the stellar mass 
was estimated as a rough mean decrease in stellar mass, over 
the entire metallicity range, between the ZAHB structures lo¬ 
cated close to the blue edge and at the center of the instabil¬ 
ity strip. The luminosity level of M evo was assumed equal to 
0.2 dex higher than the ZAHB luminosity level. Once again 
the assumed luminosity level is a rough estimate of the in¬ 
crease in luminosity typical of evolved RRLs over the entire 
metallicity range. To provide the physical structure and the 
linear eigenfunctions adopted by the hydrodynamical models, 
we computed at fixed chemical composition, stellar mass and 
luminosity level a sequence of radiative hydrostatic envelope 
models. The physical assumptions adopted to construct the 
linear models are summarized in Appendix 1, together with 
the linear blue edges based on linear hydrostatic models. The 
mean absolute bolometric magnitude and effective tempera¬ 
ture of the pulsation models approaching a stable nonlinear 
limit cycle is evaluated as an average in time over the entire 
pulsation cycle. This means that they are slightly different 


when compared with the static initial values, i.e. with the 
values the objects would have in case they were not variables. 
The difference is correlated with the luminosity amplitude and 
the shape of the light curve of the variables. The reader inter¬ 
ested in a more detailed discussion is referred to IBono et al.l 
( 1995 ). The difference in luminosity for the models located 
outside the instability strip is negligle. The difference in ef¬ 
fective temperature is negligle for models hotter than the blue 
edge and marginal for those cooler than the red edge. In Fig¬ 
ure [2] we show the location in the Hertzsprung-Russell dia¬ 
gram of a set of RRL models at fixed chemical composition 
(Z—0.0003, Y=0.245). The FU models are marked with filled 
circles, while the FOs with open circles. The black symbols 
mark pulsation models computed assuming the same stellar 
mass (0.716 M 0 ) and three different luminosity levels: the 
ZAHB (sequence A), a luminosity level 0.1 dex brighter than 
the ZAHB (sequence B) and the luminosity level of central 
He exhaustion (sequence D). The red symbols display RRL 
models computed assuming a stellar mass ~10% smaller (0.65 
M e ) than the Mzahb mass value and 0.2 dex brighter than the 
ZAHB luminosity level (sequence C). This sequence of pulsa¬ 
tion models was computed to take account of evolved RRLs. 

The adopted stellar parameters of the computed pulsation 
models are listed in Table ??. In the first three columns are 
given the metallicity, the helium content and the stellar mass 
on the ZAHB at the center of instability strip (M zahb)- The 
fourth column lists the three selected luminosity levels. For 
each adopted chemical composition, stellar mass, luminosity 
level and effective temperature, we investigated the limit cy¬ 
cle stability of RRL models both in the FU and FO mode. The 
nonlinear pulsation equations were integrated in time till the 
limit cycle stability of radial motions approached their asymp¬ 
totic behavior, thus providing robust constraints not only on 
the boundaries of the instability strip (IS), but also on the pul¬ 
sation amplitudes. 

3. NEW METAL-DEPENDENT PULSATION RELATIONS 

The correlation between pulsation and evolutionary observ¬ 
ables is rooted in several analytical relations predicting either 
the absolute magnitude (PLC relations) or the pulsation period 
(pulsation relation) as a function of stellar intrinsic parame¬ 
ters (Bono et al. 2015, in preparation). The pulsation rela¬ 
tion and its dependence on stellar mass, stellar luminosity and 
effective temperature was ca st in the form currently ad opted 
more than 40 years ago by (Ivan Albada & Bakenll971l) . To 
derive the so-called van Albada & Baker (vAB) relation they 
adopted linear, nonadiabatic, convective models. The radial 
pulsation models are envelope models, i.e. they neglect the 
innermost and hottest regions of stellar structures. This means 
that they neglect nuclear reactions taking place in the center of 
the stars and assume constant luminosity at the base of the en¬ 
velope. The base of the envelope is typically fixed in regions 
deep enough to include a significant fraction of their envelope 
mass, but shallow enough to avoid temperatures hotter than 
~ 10 7 K. This means that the computation of pulsation model 
does require the knowledge of the mass-luminosity relation of 
the stellar structures we are dealing with, van Albada & Baker 
in their seminal inve stiga tions adopted evolutionary prescrip¬ 
tions forHB stars bv llben & Rood! (119701) . 

Modern versions of the vAB relation have been derived 
using updated evolutionary models and/or nonlinear pulsa¬ 
tion models (lCoxI 119741) . More recently, they have also 
been derived using nonlinea r convective models and including 
a metallicity term (see e.g. IBono et alJ[l997l : IMarconi et alJ 
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120031: iDi Criscienzo. Marconi & Caputo 1 120041) to take ac¬ 
count of the metallicity dependence of the ZAHB luminosity 
level. In view of the current comprehensive approach in con¬ 
structing RRL pulsation models, we computed two new vAB 
relations for FU and FO pulsators. We found 

log P F = (11.347 + 0.006) + (0.860 ± 0.003) log L/L 0 + 

- (0.58 ± 0.02) log M/M q - (3.43 ± 0.01) log T e + 

+ (0.024+ 0.002) log Z (1) 

log P F0 = (11.167 + 0.002) + (0.822 + 0.004) log L/L Q + 

- (0.56 ± 0.02) log M/M q - (3.40 + 0.03) log T e + 

+ (0.013 +0.002) logZ (2) 

where the symbols have their usual meaning. The standard 
deviation for FU pulsators is 0.06 dex, while for FO pulsators 
is 0.002 dex. The decrease in the intrinsic spread for FO pul¬ 
sators was expected, since the width in temperature of the re¬ 
gion in which FOs attain a stable nonlinear limit cycle is on 
average 60% of the one for FU pulsators. 

Figure [3] shows the comparison between the computed FU 
periods and the ones predicted by the above pulsation rela¬ 
tion (black symbols) and the same comparison but adopting 
the relation bv lDi Criscienzo. Marconi & Caputo I (120041) (red 
symbols). We notice that the differences between the quoted 
pulsation relations are within the intrinsic scatter of the linear 
regressions. 

When dealing with cluster variables, one of the typical ap¬ 
proaches to improve the sample size is to fundamentalize 
the periods of FO pulsators. The classical relation adopted 
for the fundamentalization _is log / J /-v = log /Vo + 0T27 
(IDi Criscienzo. Marconi & Caputo 1 12004k iBraga et akl 12014l 
and Coppola et al. 2015, in preparation). However, this re¬ 
lation is based on a_very limited sample of observed double¬ 
mode pulsators dPetersenll 199 1) and its applicability to RRL 
has never been verified. 

Figure [4] shows the relation between FO and FU period for 
models located inside the so called OR region of the insta¬ 
bility strip: the region located between the FU blue edge and 
the FO red edge (see below), where the two modes approach 
pulsationally stable nonlinear limit cycles. The ensuing linear 
relation (black solid line) is compared with the empirical rela¬ 
tion (red solid line) obtained by using a sample of ~ 80 known 
double-mode variables identified in different stellar systems 
(Galactic globulars, dwarf spheroidals) and available in the 
literature (Coppola et al. 2015, in preparation). The classical 
relation is also shown for comparison (green solid line). Data 
plotted in this figure indicate that the new theoretical relation 
(black line) attains FU periods that are, at fixed FO period, 
slightly longer than the empirical ones (red line). On the other 
hand, the classical fix (green line) attains FU periods that are, 
at fixed FO period, slightly shorter than the empirical ones. 
We conclude that the comparison shown in Figure ^indicates 
that the classical relation is a very plausible fix over a broad 
range of metal abundances. 

4. TOPOLOGY OF THE INSTABILITY STRIP 

The current theoretical framework allows us to predict 
the approach to nonlinear limit cycle stability of the differ¬ 
ent pulsation modes. This implies the opportunity to con¬ 
strain the topology of the instability strip, i.e. the regions of 


the Hertzsprung-Russell diagram in which the radial modes 
approach a pulsationally stable nonlinear limit cycle. The 
anonymous referee noted that for the FU/FO models the ap¬ 
proach to limit cycle stability for a nonlinear system would 
imply strictly periodic oscillations. However, in the current 
theoretical framework the approach to a nonlinear limit cy¬ 
cle stability also means small changes in the mean magni¬ 
tude and effective temperature (see §2.2). An original ap¬ 
proach to compute exact periodic solutions of the nonlinear 
radiative pul sation equations was presented by Stellingwerf 
(fl97l IT981 . However, we still lack a similar relaxation 
scheme-Flouquet analysis-for radial oscillations taking ac¬ 
count of a time-dependent convective transport equation. A 
similar, but independent approach was also developed by 
iBuchler & Gounill (119841) using amplitude equation formal¬ 
ism, i.e. the temporal evolution of modal amplitudes are de¬ 
scribed by a set of ordinary differential equations. Canonical 
amplitude equations, including cubic terms, were derived to 
investigate radiative models, but nonlinear convective models 
required the inclusion of quintic terms (IBuchler et alJll999b . 
However, thei r calculation using st atic models is not triv¬ 
ial, therefore. iKo llath et al.1 (119981) . iSzabo et all (120041) and 
ISmolec & Mos kalikl (I2008al) decided to couple the solution 
of nonlinear conservation equations with amplitude equation 
formalism. The nonlinear limit cycle stability based on this 
approach is very promising, but it is very time consuming, 
since radial motions have to be analysed over many pulsation 
cycles. 

However, in the current context, we define that a radial 
mode approaches a pulsationally stable nonlinear limit cy¬ 
cle when period and amplitudes over consecutive cycles at¬ 
tain their asymptotic behavior. Stricto sensu they are not ex¬ 
actly periodic, but lato sensu they approach a_periodic be¬ 
havior. Therefore, following iBono & Stellingwerfl (119931) 
and IBono et al] d2000l) . the limit cycle in a nonlinear time- 
dependent convective regime was considered pulsationally 
stable when the differences in the pulsation properties over 
consecutive cycles become smaller than one part per thou¬ 
sands. This means that we are integrating the entire set of 
equations for a number of cycles ranging from a few hundreds 
to several hundreds. 

To support on a more quantitative basis the above defini¬ 
tion, Fig. 5 shows the nonlinear total work integral as a func¬ 
tion of integration time for three different FO (left panels) and 
FU models (right panels). They are centrally located in the 
middle of the instability strip, and constructed assuming three 
different metal abundances, stellar masses and luminosity lev¬ 
els (see labeled values). Theoretical predictions plotted in this 
figure display that the dynamical behavior, after the initial per¬ 
turbation 7 , approaches the pulsationally stable nonlinear limit 
cycle. The transient phase is at most of the order of ~ 100 
cycles (top left model), and indeed the relative changes in to¬ 
tal work after this phase are smaller than ±0.0001. To fur¬ 
ther constrain the approach to a pulsationally stable nonlinear 
limit cycle. Fig. [6] shows the bolometric amplitude as a func¬ 
tion of the integration time for the same models of Fig. 0 
Once again the luminosity amplitudes approach their asymp¬ 
totic behavior after the transition phase. After this phase the 
relative changes in bolometric amplitudes over consecutive 
cycles are smaller than ±0.001 mag. 

7 The nonlinear analysis was performed by imposing a constant velocity 
amplitude of 10 km/s both to FU and FO linear radial eigenfunctions (see 
Appendix). 
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To constrain the location of the boundaries of the instability 
strip we adopted, for each fixed chemical composition, a step 
in effective temperature of 100 K. The sampling in tempera¬ 
ture of the pulsation models becomes slightly coarser across 
the instability strip. The effective temperatures of the edges 
of the instability strip for FU and FO pulsators and for the 
adopted chemical compositions are listed in Table ??. The 
blue (red) edges are defined, for each sequence of models, 50 
K hotter (cooler) than the first (last) pulsating model in the 
specific pulsation mode. 

Figure [7] shows the boundaries for FU (solid lines) and 
FO (dashed lines) for three selected chemical compositions, 
namely Z = 0.0001 (top panel), Z = 0.001 (middle panel) 
and Z = 0.02 (bottom panel). The boundaries of the insta¬ 
bility strip, when moving from the hot to the cool region of 
the HR diagram, are the FO blue edge (FOBE), the FU blue 
edge (FBE), the FO red edge (FORE) and the FU red edge 
(FRE). The region of the instability strip located between the 
FBE and the FORE is the so called "OR region". In this re¬ 
gion the RRLs could pulsate simultaneously in the FU and 
in the FO mode. Note that we still lack a detailed knowl¬ 
edge of the physical mechanisms that drive the occurrence of 
mixed-mode pulsators (iBono et mm. This means that we 
still lack ab initio hydrodynamical calculations approaching a 
pulsationally stable nonlinear double-mode limit cycle. The 
occurrence of a narrow region of the instability strip in which 
double-mode pulsators attain a stable nonlinear limit cycle 
was suggested several years ago bv lSzabo et al.l (120041 ). How¬ 
ever, the approach adopted by these authors to deal with the 
turbulent so urce function in convectively s table regions was 
criticized by ISmolec & Moskalikl (l2008albl) . The occurrence 
of double-mode pulsators has been investig ated among clas¬ 
sical Cepheids JSmolec & Moskalikl l2008bl) and 6 Scuti/SX 
Phoenicis stars (IBono et al.l 120021) However, the modelling 
of these interesting objects is still an open problem that needs 
to be addressed on a more quantitative basis. 

It is empirically well known that mixed-mode RRLs are lo¬ 
cated in a defined region of the instability strip, between the 
long period tail of FO pulsators and the short period tail of FU 
pulsators (Coppola et al. 2015, in preparation). The above ob¬ 
servational scenario js further supporting theoretical predic¬ 
tions (iBono et al.ll 1997 ) suggesting the occurrence of mixed¬ 
mode pulsators in a narrow range in effective temperatures of 
the instability strip. 

The instability strips plotted in Figure [7] display several in¬ 
teresting features worth being discussed in more detail. 

i) - The increase in the metal content causes a shift of the in¬ 
stability strip towards cooler (redder) effective temperatures. 
This evidence was originally suggested by (Bono et alj ( Qj997l ) 
and supported by em pirical evidence (Stetson et al.l l2014t 
iFiorentino et alJ2014l) . Note that the current prediction for the 
FOREs in the metal-rich regime, need to be cautiously treated. 
There is evidence that current prediction are slightly redder 
than suggested by empirical evidence (IBono et al.lll997l) . 

ii) - The instability strip becomes, when moving from 
metal-poor to metal-rich pulsators, systematically fainter. 
Moreover, the range in luminosities covered by the three 
selected luminosity levels increases when moving from 
metal-poor to metal-rich stellar structures. The above ev- 
iden ce is a direct consequence of HB evolutionary proper¬ 
ties dPietrinferni et al.l 12004 120061). original ly brought for¬ 
ward on an empirical basis by ([Sandage 1 993). 

iii) - The region of the instability strip in which the FOs 
approach a pulsationally stable nonlinear limit cycle van¬ 


ishes at higher luminosity levels. The FOBE and the FORE 
tend to approach the same effective temperature. This point 
of the instability strip was originally called "intersection 
point” (IStellingwerfl 19751) and suggests the lack of long 
perio d FO pulsators ( Bono & Stellingwer3ll994 IBono et all 
119971). This evidence is soundly supported by observa¬ 
tions (iKunder et al.ll2013l) . In passing, we note that the cur¬ 
rent predictions and similar calculthions jwailable in the litera¬ 
ture ( Bono & Stellingwerfl 1994IBono etaJll 19971) do suggest 
the possible presence of stability isles, i.e. regions located at 
luminosities higher than the intersection point in which the 
FOs attain a pulsationally stable nonlinear limit cycle. More 
detailed theoretical and empirical investigations are required 
to constrain the plausibility of the above predictions. 

iv)- The width in temperature of the instability strip in 
which FO pulsators attain a pulsationally stable nonlinear 
limit cycle becomes systematically narrower when moving 
from metal-poor to metal-rich pulsators. The predicted trend 
appears quite clear, but we still lack firm empirical con¬ 
straints. Note that the metal-rich regime is not covered by 
cluster variables, since the most metal-rich globulars hosting 
RRLs are NGC 6388 and NGC 6441 and both of them are 
more metal-poor than Z~0.006. The empirical scenario is 
still hampered by the lack of wide spectroscopic surveys of 
Halo and Bulge RRLs. 

v')- We note that the extreme edges of the strip, namely the 
FOBE and the FRE follow a linear behaviour and this occur¬ 
rence is true for all the assumed metal contents. On the basis 
of this evidence, we derived the following analytical relations 
for the FOBE (rms = 0.003) and the FRE (rms= 0.006) as a 
function of the assumed metallicity. 

log T e F0BE = (-0.080 ± 0.008) log L/L 0 + 

- (0.012 ± 0.002) logZ + 3.957 + 0.003 (3) 


log T e FRE = (-0.084 ± 0.009) logL/L 0 + 

- (0.012 ± 0.002) logZ + 3.879 ± 0.006 (4) 

Note that the above relations suggest that the width in ef¬ 
fective temperature of the entire instability strip, among the 
different chemical compositions, is constant and of the order 
of BOOK. 

Figure [8] shows the comparison between these rela¬ 
tions and the edges_(FOBE, FRE) previously determined 
by iDi Criscienzo. Marconi & Caputo I (( 2004 ) for the three la¬ 
belled metal contents. The agreement is good and the differ¬ 
ences are within +50 K in effective temperature, that is half of 
the step in temperature adopted in the quoted grids of models. 

4.1. Pulsation properties inside the OR region 

We already mentioned in Section 4, that the modelling of 
radial pulsators that oscillate simultaneously in the FU and in 
the FO (double-mode) is still an open problem. In this con¬ 
text we define "OR region” the region of the instability strip 
in which pulsation models, after the initial perturbation, at¬ 
tain a pulsationally stable nonlinear limit cycle either in the 
FU or in the FO. This means that the approach to the nonlin¬ 
ear limit cycle does depend on the adopted initial conditions 
(linear radial eigenfunctions). From the theoretical point of 
view the models located in the OR region were adopted to 
mimic the properties of double-mode pulsators. The same 
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objects are called, from the observational point of view, RRd- 
type variables. The RRd variables play a fundamental role in 
constraining the evolutionary and the pulsation properties of 
RRLs, Indeed, the so-called Petersen dia gram (Pe tersenl 19911 
iBono et alii 19961: iBragaglia et al.l l2001). i.e. the period ratio 
between FO and FU periods ( Pfo/Pfu ) versus the FU peri¬ 
ods ( Pfu ), is a good diagnostic to constrain the pulsation mass 
and the intrinsic luminosity of RRd variables. The key advan¬ 
tage in using the Petersen diagram is that the adopted observ¬ 
ables are independent of uncertainties affecting the distance 
and the reddening correction of individual objects. Moreover 
and even more importantly, theory and observations indicate 
that both the mean magnitudes and colors of the two modes 
are within the errors the same (ISoszynski et al.|[2014l) . The 
above evidence and the pulsation relations discussed in §3 in¬ 
dicate that the Petersen diagram can be soundly adopted to 
constrain the actual mass of RRd variables. 

A detailed theoretical investigation of double-mode RR 
Lyrae in the Petersen diagram was provided by lPopielski et al.l 
(l2000l) . They investigated in detail the dependence of the pe¬ 
riod ratio on intrinsic parameters (stellar mass, luminosity, ef¬ 
fective temperature, metal abundance). Moreover, they also 
performed a detailed comparison with RRd in several globular 
clusters and nearby dwarf galaxies. In particular, they found 
that Large Magellanic Cloud double mode variables cover a 
modest range in metallicity (-1.7<[Fe/H]<-1.3). Moreover, 
the spread in period ratio at fixed iron abundance could be ex- 
plained as a spread in stellar mass. Theoretic al (IBono et al.l 


[1996tjKovacs & Walker 

1999; KovacsUjOOO) and empiri- 

cal (Beaulieu et al. 1997 

Alcock et al. 1999; Soszvnski et al. 


1201 11120141) evidence indicates that more massive RRd vari¬ 
ables attain, at fixed FU period, larger period ratios. 


Moreover, theory and observations suggest that metal-rich 
RRd variables show, when compared with metal-poor objects 
shorter FU periods and smaller period ratios (|Bragagl ia et al.l 

l2001tIPi Criscienzo et al .11201 mSoszvhski etal .112014l)~ 

To take account of both the dependence on the stellar mass 
and on the metallicity, we computed the period ratios of the 
pulsation models located in the ”OR region” for each assumed 
chemical composition. Figure [9] shows the distribution of the 
quoted models in a 3D logarithmic plot including the period 
ratio ( Pfo/Pfu X the stellar mass in solar units and the metal 
content. Data plotted in this figure display a well defined cor¬ 
relation. Therefore, we performed a Least Squares linear re¬ 
gression and we found the following analytical relation: 

logM/Mo = -0.85(±0.05) - 2.8(±0.3 )log(P FO /P FU ) + 

-0.097(+0.003 )logZ (5) 

where the symbols have their usual meaning. The above 
relation, with a rms of 0.004, provides the unique opportunity 
to constrain the pulsation mass of double-mode pulsators on 
the basis of their period ratios and metal-contents (Coppola et 
al. 2015, in preparation), observables that are independent of 
uncertainties affecting individual distances and reddenings. 

5. THE PERIOD-RADIUS RELATION 

Recen t improvemen ts in inte r ferom etric measure¬ 
ments dKervellal 2008 : lErtel et akl 1 2014 1) provided the 
opportunity to measure the diameter of several evolved radial 
variables (classical Cepheids, Kervella et al. 2001; Mira, 
Milan-Gabet 2005). The next generation of optical and 
NIR interferometers (Nordg ren et al.l l2000t iGallenne et al.l 


120141 ) will allow us to measure the diameter of nearby field 
RRLs. Moreover, recent advancements in the application 
of the Infrared Surface Brightness (IRSB) method to RRLs 
is providing new and accurate measurements of RRL mean 
radii. 

Therefore, we decided to use the current sequences of non¬ 
linear, convective models to constrain the Period-Radius (PR) 
relation of FU and FO RRLs. Figure [TO] shows in a loga¬ 
rithmic plane the mean radius of FU pulsators over the entire 
set of chemical compositions versus the FU period. We per¬ 
formed a linear Least Squares regression over the entire set of 
models and we found the following PR relation for FU pul¬ 
sators (black line): 

logR/R G = 0.866(+0.003) + 0.55(+0.02) logP (6) 

with a standard deviation of 0.03 dex. The red and the green 
solid lines display the predicted PR relation bv lMarconi et al.l 
(120051) . based on similar RRL models and the extrapolation 
to shorter periods of the predicted PR relation of BL Her- 
culis variables, provided by [Marconi & Di Criscienzd l d2007l) 
(MDC07), using the same theoretical framework adopted in 
this investigation. The standard deviations of the individual 
PR relations plotted in the top left corner of the same figure 
indicate a good agreement among the different predicted PR 
relations over the entire period range. 

The difference between the new PR relation and the one de¬ 
rived by dMarconi et al.ll2005l) is the consequence of the dif¬ 
ferent assumptions concerning stellar masses and luminosities 
adopted for the different chemical compositions. In the cur¬ 
rent approach, we adopted evolutionary prescriptions, while 
in ( Marconi et al.l2005l) the grid of pulsation models was con¬ 
structed by adopting, for each chemical composition and stel¬ 
lar mass, a fixed step in luminosity level. In spite of the 
different approach adopted in selecting the evolutionary pa¬ 
rameters, the agreement is quite good. The difference be¬ 
comes slightly larger only in the BL Herculis regime, i.e. for 
log P > 0. 

Note that the period range adopted in the above figure is 
larger than the typical range of RRL stars (~ 0.2 < P < 1.0 
days, Marconi et al. 2011). To validate the above theoret¬ 
ical scenario Figure [TO] also shows the comparison with the 
empirical PR relation provided by iBurki & Mevlanl (119861) 
(dashed green line) for both RRL and BL Herculis variables. 
The agreement is once again quite good over the entire period 
range. 

The similarity of both predicted and empirical PR rela¬ 
tions for RRL and BL Herculis supports earlier suggestions 
concerning the tight evolutionary and pulsation correlation 
of these two classes of evolved low-mass radial variables. 
This applies not only to the PR rela tions (IBurki & Mevlanl 


1986 

Marconi & Di Criscienzo 2003), but also to the PL rela- 

tion 

Caputo et al. 2004; Matsunaga et al.][2009; Ripepi et al.l 

2015 

). In this context, the BL Herculis are just the evolved 


component of RRL stars, i.e. HB stellar structures that in their 
off-ZAHB evolution cross the instability_strip at higher lumi¬ 
nosity levels than typical RRLs dMarconi et ali201 1 ). In spite 


of this indisputable similarity between RRL and BL Herculis 
a word of caution is required. According to the above evo¬ 
lutionary framework the BL Herculis evolve from the hot to 
the cool side of the instability stripe. This means that they 
require a good sample of hot HB stars. However the HB mor¬ 
phology, i.e. the distribution of HB stars along the ZAHB, 
does depend on the metallicity. The larger is the metallic- 
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ity, the redder becomes the HB morphology ( Cast ellanllll983b 
lRenzinilll983l) . This means that the probability to produce BL 
Herculis becomes vanishing in metal-rich stellar systems. On 
the other hand, we have evidence of RRL stars with iron abun¬ 
dances tha t are either solar or even super-solar in the Galac¬ 
tic Bulge (I Walker & Terndrupll 19911) . This would imply that 
the above similarity between RRL and BL Herculis variables 
might not be extended over the entire metallicity range. 

Predicted mean radii plotted in Figure [TO] display, at fixed 
period, a large intrinsic dispersion. To further constrain this 
effect and to estimate the sensitivity of the PR relation on 
the metal content, we also computed new Period-Radius- 
Metallicity (PRZ) relations for both FU and FO pulsators. We 
found 

log R/R 0 = 0.749(+0.006) + 0.52(±0.03) log P + 

- 0.0391+0.006) logZ (7) 

with a standard deviation of 0.006, for FU models and 

log R/R 0 = 0.87(±0.03) + 0.601+0.06) log P + 

- 0.033(±0.005) logZ (8) 

with a standard deviation of 0.006, for FU models and 0.003 
forFO models. 

The models and the new PR relations are plotted in the top 
(FU) and in the bottom (FO) panel of FigureQj] Again, simi¬ 
lar relations by iMarconi et al.l 1 2005 ) are shown for compari¬ 
son. As expected the PRZ relations display more tight corre¬ 
lations than the classical PR relation. Thus suggesting a clear 
dependence on the metal content. The above evidence indi¬ 
cates that the PRZ relation is a powerful tool to constrain in¬ 
dividual radius estimates for RRL of known period and metal¬ 
licity with a precision of the order of 1%. 

6. PREDICTED LIGHT CURVES 

One of the key advantages in using nonlinear, convective 
hydrodynamical models of RRL stars is the possibility to pre¬ 
dict the variation of the leading observables along the pulsa¬ 
tion cycle. Among them the bolometric light curves play a 
fundamental role, since mean magnitudes and colors do de¬ 
pend on their precision. In the current investigation, we com¬ 
puted limit cycle stability for both FU and FO RRLs cover¬ 
ing a broad range in chemical compositions and intrinsic pa¬ 
rameters (stellar mass, luminosity levels). The different se¬ 
quences of models do provide an atlas of bolometric light 
curves that can be adopted to constrain not only the mean 
pulsation properties of RRLs, but also the occurrence of sec¬ 
ondary features (bumps, dips) along the pulsation cycle. The 
Figures[T2]and[T3]show predicted FU and FO bolometric light 
curves for the sequence of models at Z = 0.0006, Y = 0.245, 
M = 0.67M 0 , logL/L G = 1.69 s . 

6.1. Transformation into the observational plane 

The bolometric light curves discussed above have been 
transformed into the most popular optical (UBVRI) and NIR 
(JHK) bands using static mode l atmospheres dBono et al.1 
fl99llCastelli. Grafton & Kuruczlll997alhll . 

Once the bolometric light curves have been trans¬ 
formed into optical and NIR bands, the magnitude- 

8 The bolometric light curves for the other chemical compositions are 
available upon request to the authors. 


averaged and intensity-averaged mean magnitudes and col¬ 
ors can be derived for the entire set of stable pulsa¬ 
tion models. The intensity-averaged magnitudes are the 
most used in the literature, since they overcome differ¬ 
ent thorny problems connected with the shape of the light 
curves dBono & Stellingwerf)ll993i) 9 . They are presented in 
the following. 

6.2. Mean magnitudes and colors 

The Tables ?? and ?? give the intensity-averaged mean 
magnitudes for the entire grid of FU and FO models. For 
each set of models the Tables list the chemical composition 
(metal and helium content) together with the intrinsic pa¬ 
rameters (stellar masses and luminosity levels) adopted for 
the different sequences. The first three columns list, for 
every pulsation model, the effective temperature, the loga¬ 
rithmic luminosity level and the logarithmic period. The 
subsequent eight columns give the predicted Johnson-Kron- 
Cousins (UBVRI) and NIR (JHK, 2MASS photometric sys¬ 
tem) intensity-averaged mean magnitudes. 

When moving the instability strip from the HR diagram into 
the observational planes, the predicted FOBE and FRE colors 
can be correlated with the absolute magnitude and the metal 
content according to the relations given in Table ??. We note 
that in the case of K — (J — K) the predicted boundaries do not 
depend on metallicity. 

7. METAL-DEPENDENT PERIOD-LUMINOSITY AND 
PERIOD-WESENHEIT RELATIONS 

The current grid of nonlinear, convective RRL models pro¬ 
vide the unique opportunity to build a new and independent 
set of metal-dependent diagnostics to determine RRL dis¬ 
tances. We focussed our attention on metal-dependent PL 
(PLZ) relations and on the metal-dependent PW relations 
(PWZ). This approach was adopted for the entire set of op¬ 
tical, optical-NIR and NIR intensity-averaged mean magni¬ 
tudes. 

Note that to derive the PLZ and the PWZ relations for FU 
pulsators we included, for each assumed chemical composi¬ 
tion, sequences A,B and C of models (see Table 2 and Figure 
2). We neglected the sequence D models, since these lumi¬ 
nosity levels are too bright for the bulk of RRLs. They are 
more typical of radial variables at the transition between BL 
Herculis and W Virginis variables. 

To derive the PLZ and the PWZ relations for FO pulsators 
we included the two faintest luminosity levels (sequences A 
and B). We neglected the brighter sequences C and D, since 
they are typically brighter than the ’’Intersection Point” (see 
§4). 

7.1. The PLZ relations 

The predicted PLZ relations were derived for FU and FO 
pulsators covering the entire set of chemical compositions 
(Z=0.0001-0.02). We performed several linear Least Squares 
relations of the form: 

M x = a + blogP + c[Fe/H] (9) 

Together with the independent PLZ relations for FU and 
FO pulsators we also estimated an independent set of’’global” 
PLZ relations using simultaneously FU and FO pulsators. The 

9 The magnitude-averaged quantities are available upon request to the au- 
thors. 
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periods of the latter group were fundamentalized using the 
classical relation (see also §3). This is the typical approach 
adopted to improve the sample size, and in turn, the precision 
of cluster/galaxy distances based on RRLs (iBraga et al.ll20l4 
Coppola et al. 2015, in preparation). The classical relation is 
log Ppu = log Pk) +0.127 and the symbols have their usual 
meaning. 

Theo retical (IBono et al. l2001|JCatglungtal.l l2004l) and em¬ 
pirical dBenko etalJ 1201 1 : IBraga et alJl2014l) evidence indi¬ 
cates that RRLs do not obey to a well defined PL relation 
in the blue ( UBV ) bands. The slope of the PL relation be¬ 
comes more and more positive when moving from the R to 
the NIR bands. The slope is the consequence of a signif¬ 
icant variation in the bolometric corrections when moving 
from the blue (short periods) to the red (long periods) edge 
of the instability strip. This change is minimal in the blue 
bands and becomes larger than one magnitude in the NIR 
bands (iBono et al.ll2003l) . These are the reasons why we de¬ 
rived new PLZ relations, only for the RIJHK bands. The 
coefficients of the RIJHK PLZ relations are listed in Table ?? 
together with their uncertainties and standard deviations. 

The left panels of Figure [14] show the PLZ relations for FU 
and FO pulsators, while the right ones the global PLZ rela¬ 
tions. Lines of different color display relations with metal 
contents ranging from Z = 0.0001 (black) to Z = 0.02 (pur¬ 
ple). The current models soundly support the empirical evi¬ 
dence dBenko et alJl201 it ICoppola et al.ll20lTb that the slope 
of the PL relations becomes steeper when moving from the 
optical to the NIR bands. Indeed, the slope increases from 
1.39 in the R-band to 2.27 in the K -band. Moreover, the stan¬ 
dard deviation of the PLZ relations decreases by a factor of 
three when moving from the optical to the NIR bands. This 
means that the precision of individual distances, at fixed pho¬ 
tometric error, increases in the latter regime. The above trends 
are mainly caused by the temperature dependence of the bolo¬ 
metric correction as a function of the wavelength. 

Moreover, optical and NIR absolute magnitudes of metal- 
poor RRLs are, at fixed period, systematically brighter than 
metal-rich ones. The difference is mainly due to evolutionary 
effects. The ZAHB luminosity at the effective temperature 
typical of RRLs (log 7’ e =3.85) decreases for an increase in 
metal content dPietrinferni et al.ll2013l) . Note that the metal- 
licity dependence is, within the errors, similar when moving 
from the optical to the NIR bands. Indeed the coefficient of 
the metallicity term attains values of the order of 0.14-0.18 
dex. 

Figure[l5]shows the comparison between predicted and em¬ 
pirical K -band PL relations available in the literature. The top 
panel shows the comparison for the most metal-poor chemi¬ 
cal composition (Z = 0.0001), while the bottom panel for a 
metal-intermediate chemical composition (Z = 0.001). Note 
that in the comparison we only included PL relations tak- 
ing account either of a metallicity term (IBono et al.l 1200 It 
Catelan et al.ll2004h or of the HB morphol ogy dCassisi et all 


20041) or based on an empirical approach (IBono et al.l 20031: 


Sollima et alj|2006b . 


It is worth mentioning that the PLZ relation provided 
by IBono et ah] ( 2003 ) is based on field and cluster RRLs for 
which the individual distances were estimated by using the 
NIR surface-brightness method (IStorm et al.l 1 1994 ft. On the 
other hand, the PLZ relation provided bv ISollima et al] (120061) 
was estimated using RRLs from 16 calibrating GCs to fix 
the slope of the relation and RR Lyr itself to fix the zero- 
point of the relation. The above comparison indicates that 


the current K -band PL relations_agree quite well with sim¬ 
ilar predicted ( [Catelan et al.l 120041) (blue lines) and empiri¬ 
cal (ISollima et al.l 120061) (cyan lines) PL relations. The dif¬ 
ference is on average smaller than 1 cr both for metal-poor and 
metal-intermediate chemical compositions. 

The current predictions attain intermediate sl opes w hen 
compared with similar prediction s bv IBono et al.l (120011) (red 
lines) and by IBono et aD (120031) (green lines). The differ¬ 
ences are caused either by different assumptions concern¬ 
ing the ML relations_adopted to construct the grid of pul¬ 
sation models ( Bono et all 120011) or by different assump¬ 
tions concerning the absolute calibrators and/or the metallic¬ 
ity scale (IBono et al.ll2003l) . 

The above discussion is mainly focussed on the compar¬ 
ison between predicted and empirical slopes of the K -band 
PL relations. A glance at the PL relations plotted in top 
panel of Figure [T5Tindicates that the zero-point of the current 
metal-poor prediction is ~0.1 mag brighter than the empir¬ 
ical PL relation by ISollima et al.l 020061) . However, the dif¬ 
ference becomes marginal in the metal-intermediate regime 
(bottom panel). The difference might be the consequence 
of the adopted zero-point. The empirical K -band PL rela¬ 
tion is rooted on RR Lyr itself, which is a metal-intermediate 
([Fe/H]=-1.50+0.13 dex; Braga et al. 2014) field RRL. 

7.2. Dual band PWZ relations 

The PLZ relations have several advantages in estimat¬ 
ing individual distances. However, they are prone to un¬ 
certainties affecting reddening corrections. This problem 
is strongly miti gated when moving into the NIR and MIR 
regime (IMadore & Freedman! 120121) . but still present. The 
problem becomes even more severe in regions affected by dif¬ 
ferential reddening. T o overcome this prob l em it was sug- 
gested ~40 years ago (Ivan den Berghl 119751 : iMadorel 119821) 
to use reddening free pseudo-magnitudes called ’’Wesenheit 
magnitudes”. They are defined as: 

W(X, Y) - M x + %(M X - M y ) (10) 

where the coefficient of the color term - £ - is the ratio be¬ 
tween the selective absorption in the X band and the color ex¬ 
cess in the adopted color. It is fixed according to an assumed 
reddening law that in our case is the lCardelli et al.l d 1989b law. 
The coefficients of the color terms are listed in the third col¬ 
umn of Table ??. 


Pros and cons of both optical and NIR PW relations have 
been widely discussed in the literature (llnno et all 12013 ). 
However, they have been typically focussed on classical 


Cepheids (Caputoetal. 2000; Fiorentino et al. 

20021 [2007; 

Marconi et al. 

2005; Ngeow & Kanbur 

2005 

; Bono et al. 

2010: Storm et alJ2011; Ripepi et al. 2012 

). We briefly sum- 


marise the key issues connected with optical and NIR PW 
re lations. The read er is referred to Braga et all (120141) and 
to llnno et al.l (120131) for a more detailed discussion. The pros 
are the following: a) Individual distances are independent of 
reddening uncertainties, b) The PW relation mimic a PLC 
relation, since they include a color term. This means that in¬ 
dividual distances are more precise when compared with dis¬ 
tances based on PL relations, c) The metallicity dependence 
for classical Cepheids app ears to be vanishing for optical- 
NIR and NIR PW relations (ICaputo et alJ2000tlMarconi et alJ 
120051: llnno et alJl2013l) . The main cons are the following: a) 
The PW relation assume that the reddening law is universal. 
However, the use of bands with a large difference in central 
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wavelength mitigates the problem (llnno et al .112013l) . b) Two 
accurate mean magnitudes are required to provide the dis¬ 
tance. 

To investigate the properties of optical and NIR PW in 
the RRL regime we computed new FU, FO and global PW 
relations for the different combinations of the above five 
bands. They were derived using the entire metallicity range 
(Z=0.0001-0.02). The coefficients, their errors and standard 
deviations are listed in Table ??. The panels of Figures [16] 
shows from top to bottom the PW relations based on optical 
bands. The color coding of the different lines is the same as 
in Figure [14] A glance at the predictions plotted in this fig¬ 
ure discloses an interesting result. The PW(V.B-V) relations 
plotted in the top panels show a minimal, if any, dependence 
on the metal content. The coefficients of the metallicity term 
for FU, FO and global PWZ relations are, within the current 
uncertainties, similar and vanishing (Table ??). This finding 
is at odds with similar PWZ relations for classical Cepheids. 
The metallicity dependence for these objects is larger in the 
optical regime and becomes smaller either in the NIR or in 
the optical-NIR regime. The difference is mainly due to the 
higher range in surface gravities and effective temperatures 
covered by RRLs when compared with classical Cepheids. It 
is worth mentioning that the above finding is even more ap¬ 
pealing if we account for the fact that RRLs cover a range 
in metallicity that is at least one dex wider than the range 
covered by classical Cepheids ([Fe/H]^ -2.5 - +0.30 vs 
[Fe/H]^ -1.5 - +0.30; Pedicelli et al. 2009; Romaniello 
et al. 2008; Dambis et al. 2014). 

The above findings bring forward two important conse¬ 
quences: a) The PW(V.B-V) relation are robust distance in¬ 
dicators for RRL sharing simultaneously the advantages to be 
reddening free and almost independent of the metal content, 
b) The use of the PW(V.B-V) relation together with indepen¬ 
dent optical and NIR magnitudes can also provide tight con¬ 
straints on individual RRL reddenings. The main cons is, once 
again, the dependence on the assumed reddening law. 

To constrain this effect we adopted two independent 
redd ening laws, namely the ones by iMcCalll (120041) 
and IFitzpatrick & Mas~sal (120091) . We found that the 
color coefficient of the PW(V, B-V) relation change from - 
3.07 (IMcCalll l200l to -3.36 (IFitzpatrick & Massal l2009h . 
This difference causes a difference in distance of the order 
of 0.2%. 

The PW relations plotted in Figures [16] and listed in Ta¬ 
ble ?? indicate that the coefficient of the metallicity term be¬ 
comes of the order of 0.1 in the other optical bands, while 
it increases to 0.15-0.20 in the optical-NIR and in the NIR 
bands. Thus suggesting a similar dependence of the optical- 
NIR and NIR regimes on the metal content. 

However, the standard deviation has a different behavior. It 
is on average a factor of two smaller in the NIR-bands that in 
the optical ones. The smallest values are attained by optical- 
NIR PW relations in which they are of the order of a few hun¬ 
dredths of magnitude. The above evidence indicate that ul¬ 
timate precision of individual distances is a balance between 
photometric and spectroscopic uncertainties. Thus suggest¬ 
ing that the error budget is the main criterium in the selection 
of the appropriate PWZ relation to estimate individual RRL 
distances. 

To validate the above PWZ relations we compared the cur¬ 
rent predictions with similar mass-dependent PW relations 
available in the literature. In particular, we focussed our at¬ 
tention on the BV and the VI filter combinations provided 


by iDi Criscienzo, Marconi & Caputo I (120041) . Note that the 
latter predictions rely on models covering a narrower metal¬ 
licity range (Z < 0.006). Figure [19] shows the compar¬ 
ison for the BV filters in two different metal-abundances. 
For each selected metal abundance we adopted the stellar 
mass predicted by HB evolutionary models for a ZAHB 
structure located at log T e = 3.85 (see Table ??). These 
mass values and the the mass-dependent relations provided 
by Pi Criscienzo. Marconi & Caputo I (120041) were adopted 
for the two labelled metallicities. The agreement is excellent 
for Z = 0.001 (right panel) and within lcr for the more metal- 
poor, Z = 0.0001 (left panel), chemical composition. The 
difference is mainly caused by a difference in the zero-point, 
and in turn, in the adopted ML relation. The slopes attain both 
in the metal-poor and in the metal-rich regime similar values. 

The comparison of the VI filters is of little use, 
since the coefficient of the color term adopted 
by IDi Criscienzo. Marconi & Caputo I (120041) is slightly 
different than the current one (1.433 vs 1.54) 10 . To overcome 
the problem we adopted their coefficient and we found that 
both the zero-points and slopes agree quite well in the 
metal-poor and in the metal-rich regime. 

7.3. Triple band PWZ relations 

The dual band PW relations have been widely used in the 
recent literature in dealing with RRL and Cepheid individual 
distances. During the last few years it has been suggested 
by iRiess et ail ( 2011 ). on an empirical basis, to use triple 
bands PW relations. These new diagnostics can be defined 
as 

W(X , y, W) = Mx + t(My - M w ) (11) 

where the coefficient of the color term -ij— is the ratio be¬ 
tween the selective absorption in the X band and the color 
excess in the adopted (T - W) color. In the quoted paper 
the authors adopted the F555W and the L814 W bands plus 
the NIR TT60W band in the WFC3 HST photometric system 
and they applied the new PWZ relation applied to Classical 
Cepheids in external galaxies. The same approach, but based 
on nonlinear, convective Cepheid models was also adopted 
bv iFiorentino et all (120131) . More recently, triple bands PW 
relations have been derived by Braga et al. (2014) for optical, 
optical-NIR and NIR PW relations of RRLs in M4. 

According to the quoted authors, the key advantages in us¬ 
ing triple bands PW relations to estimate individual distances 
are the following: i) they have smaller intrinsic dispersions 
when compared with dual band PW relations; ii) they are less 
prone to possible systematics in using the same magnitude in 
the color term. The main drawback is the need of accurate 
mean magnitudes in three different bands. 

Here we derive for the first time triple band PWZ relations 
using RRL pulsation models. They are listed in Table ?? and 
shown in Figure l20l for different combinations. Note that we 
focussed our attention on triple PWZ relations including a 
NIR magnitude and an optical color. We found that the co¬ 
efficients of the color term are typically smaller than 0.7 and 
even smaller than 0.2 for colors including bands with large 

10 The difference is caused by a different assumption concerning the cen- 
tral wavelength of the Johnson-Kron-Cousins /-band adop t ed to e stimate the 
selective absorption from the reddening law (Cardelli et alJU989h . In partic¬ 
ular, IDi Criscienzo, Marconi & Caputo I (20041) adopted A/ = 7940 A, while 
we adopted Aj = 8060 A (Besse TFim 
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differences in central wavelengths (B - R, B - I). The coeffi¬ 
cients of the metallicity term attain on average values similar 
to the dual bands PWZ relations. The standard deviations of 
the quoted triple bands PWZ relations are, once again, similar 
to the standard deviations of dual bands PWZ relations. 

7.4. Metal-independent PWZ relations 

We found that a few dual and triple band PWZ relations 
listed in Tables ?? and ??, have coefficients of the metal¬ 
licity term that are smaller than 0.05 mag/dex. This im¬ 
plies a vanishing dependence on the metal content, since a 
difference of 1 dex in metal content would imply a differ¬ 
ence in distance modulus at most of the order of 0.05 mag. 
Therefore, we decided to compute for the same PWZ rela¬ 
tions and new set of PW relations that neglect the metallic¬ 
ity dependence. The zero-points and the slopes of the metal- 
independent relations are listed in Table ??. The new PW re¬ 
lations show standard deviations that are only slightly larger 
than the metal-dependent PWZ relation. Thus further sup¬ 
porting the marginal role played by the metal content over the 
entire period range. This theoretical evidence, once validated 
on an empirical basis, opens the path to a new diagnostic to 
estimate individual distances of RR Lyrae for which the metal 
content is not available. 

7.5. Uncertainties affecting the coefficients ofPLZ and PWZ 

relations 

The new theoretical framework for RR Lyrae stars we 
are developing does depend on the physical assumptions 
adopted in the treatment of turbulent convection. We use a 
nonlinear, nonlocal, time-dependent approach to deal with 
convective transport. However this treatment, and similar 
approaches available in the literature (iSmolec & Moskahkj 
l2008ah ISmolec & Mosl<al ikll201 Ol) . do rely on a free param¬ 
eter, the so-called mixing length parameter. In the current 
approach, we adopted a mixing length parameter equal to 

1.5. Plausible changes of this parameter mainly affect the 
pulsation amplitudes and to a minor extent the boundaries 
of the instability strip. Detailed calculations concerning 
the dependence of nonlinear observables on the efficiency 
of the convective transport have already been discussed by 
iDi Criscienzo, Marconi & Canuto I d2004l) . The treatment 
of the turbulent convection also relies on the use of a few 
other free parameters. They have been fixed according 
to the value of the mixing length param eter following 
the prescriptions provided b\ lBono & Stellingwerfl (119941) . 
The adopted values have been validated fitting the light 
curves of field and cluster RR Lyrae (IBono et al . 2000; 
Marconi & Clementinil 120051; [M arconi & DeglTnnocenti 
20071) . Bump Cepheids (Bono et all l2002cl) . the prototype 
6 Cephei dNatale et all 1200 81) and cl assical Cepheids in 
eclipsing binaries (iMarconi et al.l 1201 3t) . Moreover, we do 
allow the convective flux to attain negative value s at the 
boundaries of convective stable regions ( Bono et al, 2000 )■ A 
similar approach was also adopted by (ISmolec & Moskahkl 
l2008ah . 

Bolometric light curves are transformed into the ob¬ 
servational plane using bolometric corrections and color- 
temperature (CT) relations predicted by static stellar atmo¬ 
sphere models. This means that we are assuming a quasi¬ 
static approximation in transforming predicted observables 
(static vs effective surface gravity). A proper treatment does 
requ ire the detailed solution o f the radiative transfer equa¬ 
tion ( iDorfi & Feuchtingeill999l) . In spite of the limitations of 


the current approach, we found that transformations based on 
an independent set of stellar atmosphere models (PHOENIX 
iKucinskas et al.ll2006l) do provide very similar results. In this 
context, it is worth mentioning that a detailed comparisons 
between predicted and observed mean effective temperatures 
of RR Lyrae stars indicate a difference of the order of 150 K 
(ICacciari et alj|2000h . Moreover, we still lack a detailed com¬ 
parison between color-temperature relations based on static 
and hydrodynamic al atmosphere models. 

In the current investigation we adopted, following 
iPietrinferni et al.l (12006). a (/-enhance d chemical mixture . 
Moreover, we are adopting, following Cas sisi et all (120031) . 
a primordial helium abundance of Yp- 0.245, consistent 
with measurements of the cosmic microwave backg round 
dPrvke et al .1120021 lAde et ai1l2014t iHinshaw et alJl2013l) and 

a helium-to-metal enrichment ratio of 1.4. Larger helium 
abundances, at fixe d metal content, affect the pulsation pr op- 
erties of RR Lyrae (IBono et alJI 1995t IMarconi et al.ll201 ll) . A 
systematic investigation of the impact that helium abundance 
has on pulsation observables will be investigated in a forth¬ 
coming paper. 

Finally, the evolutionary predictions adopted in this inves¬ 
tigation are also affected by uncert ainties in the adopted in¬ 
put physics (iCassisi et al.ll 19981 12007 ) and on the treatment of 
mass loss, atomic diffusion, extra-mixing and neutrino losses. 
The reader interested in a more detailed analysis of the uncer- 
tai nties affecting the pr e dicted ZAHB l u minosity is r eferre d 
t o ICassisi et ITT(fT998h . LSalarisI (f2()Tl . lYalle et all (f2()Tl . 
lVandenBergl(l201 31) . 

8. SUMMARY AND FUTURE REMARKS 

We present a comprehensive theoretical investigation of the 
pulsation properties of RRL stars. To provide a homoge¬ 
neous and detailed theoretical framework to be compared with 
the huge photometric and spectroscopic data sets that are be¬ 
coming available in the literature we compute a large grid of 
nonlinear, convective hydrodynamic al models of RRL stars. 
The RRL models were constructed assuming a broad range of 
metal abundances (Z = 0.0001-0.02) and at fixed helium- 
to-metals enrichment ratio (AT/AZ=1.4). As a whole we 
computed »420 nonlinear hydrodynamical models. Among 
them »300 display a pulsationally stable nonlinear limit cy¬ 
cle, while ~60 experience a mode switching. The latter group 
approaches, after a transient phase, a nonlinear limit cycle that 
is different from the initial perturbed linear radial eigenfunc¬ 
tion. Moreover, «60 models quench radial oscillations be¬ 
cause they are located outside the instability strip. They are 
either hotter than the FOBE or cooler than the FRE. The main 
difference of the current grid when compared with similar cal¬ 
culations available in the literature is that for each fixed chem¬ 
ical composition, the stellar mass and the luminosity levels 
adopted to construct pulsation models were fixed according to 
detailed central He burning HB evolutionary models. In par¬ 
ticular, for each fixed chemical composition we adopted two 
different stellar masses to take account of RRL located either 
in the proximity of the ZAHB (sequence A in Figure 2) or 
crossing the instability strip, from hot to cool effective tem¬ 
peratures, at higher luminosity levels (sequence C in Figure 
2). Indeed the former models were constructed by assuming 
three different luminosity levels (A, B, D) to take account of 
the off-ZAHB evolution till central He exhaustion. 

To provide a thorough analysis of the topology of the RRL 
instability strip, as a function of the metal content we com¬ 
puted the pulsation stability for both FU and FO pulsators. 
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The calculations were extended in time till the individual 
models approached limit cycle stability and we could con¬ 
strain their modal stability. The main results of the above 
theoretical framework are the following. 

• Pulsation properties. The current theoretical framework 
allowed us to provide detailed predictions of a broad range of 
observables (amplitudes: luminosity, velocity radius, effective 
temperature; mean magnitudes, velocities, radii). In particu¬ 
lar, we investigated their dependence on the metal content. 

• Modal stability. We provided a detailed mapping of the 
RRL instability strip as a function of the metal content. We 
found that an increase in metal content causes a systematic 
shift of the instability strip towards redder (cooler) colors. 
This confirms previous findings by our group. 

• Pulsation relations. We provided accurate pulsation re¬ 
lations (van Albada & Baker relations) for both FU and FO 
pulsators. The key advantage of the current relations is that 
they rely on a homogenous evolutionary and pulsation frame¬ 
work. 

• Pulsation masses of double mode pulsators. The topology 
of the instability strip allowed us to constrain the so-called 
”OR” region, i.e. the region in which both FU and FO pul¬ 
sators attain a pulsationally stable nonlinear limit cycle. Mod¬ 
els located in this region were adopted to mimic the properties 
of double mode pulsators. We derived a new analytical rela¬ 
tions to constrain the mass of double mode pulsators using 
period ratios and metal content. 

• Period radius relations. We derived new Period-Radius- 
Metallicity (PRZ) relations for FU and FO pulsators. They 
agree quite with previous PR and PRZ relations. 

• Transformations into the obserx’ational plane. Bolo- 
metric magnitudes and effective temperatures were trans¬ 
formed into the observational plane using bolometric 
corrections and color-temperature relations provided 
by Castelli. Gratton & Kuruczl (11997allbh . This means that 
we provide intensity-averaged mean magnitudes and colors 
together with luminosity amplitudes in the most popular 
optical and NIR bands ( UBVRIJHK ). 

• RR Lyrae as distance indicators. Homogeneous predic¬ 
tions concerning optical/NIR mean magnitudes allowed us to 
compute new diagnostics to estimate individual distances of 
Galactic and Local Group RRL stars. 

Period-Luminosity-Metallicity relations --We 
derived new accurate PLZ relations for RRL stars in opti¬ 
cal/NIR bands (RUHR). We confirm that RRL stars do not 
obey to PLZ relations in the blue regime. The lack of a well 
defined slope is the minimal dependence of the bolometric 
correction, in these bands, on the effective temperature. 

11 http://www.eso.org/public/teles-instr/e-elt.html 

12 http://www.tmt.org/ 

13 http://www.gmto.org/ 


Period-Wesenheit-Metallicity relations -- We 
derived new accurate PWZ relations for RRL stars in opti¬ 
cal, optical-NIR and NIR bands. The key advantages of the 
PWZ relations is that they are independent of reddening un¬ 
certainties, moreover, they also show smaller intrinsic disper¬ 
sion when compared with similar PLZ relations. The latter 
feature is the consequence of the inclusion of a color term 
mimicking a PLC relation. The main drawback is that they 
rely on the assumption that the reddening law is universal. 
However, theoretical and empirical evidence indicate that op¬ 
tical/NIR PWZ relations are less prone to secondary features 
of the reddening law. To fully exploit the use of the PWZ rela¬ 
tions as distance indicators we computed both dual and triple 
band PWZ relations. The latter appear very promising, but 
they do require three independent mean magnitudes. Finally, 
we found that the predicted PW(V.B-V) relations are almost 
independent of the metal content. 

The current investigation is the first step of a large project 
aimed at constraining the pulsation properties of RRL as a 
function of chemical composition and ages of the progenitors. 
We plan to investigate the dependence on the helium con¬ 
tent in a forthcoming investigation. Moreover, we also plan 
to transform the current predictions from the UV to the MIR 
using homogeneous sets of stellar atmosphere models. This 
new theoretical scenario shall pave the road for a massive use 
of RRL stars as distance indicators. The above plan appears 
even more compelling in waiting for the first Gaia data release 
together with the advent of new space (JWST) and ground- 
based observing facilities. In this context the Extremely Large 
Telescopes (European-ELT [E-ELT] 11 , the Thirty Meter Tele¬ 
scope [TMT] 12 and the Giant Magellan Telescope [GMT] 13 ) 
will play a crucial role, since they are going to resolve indi¬ 
vidual HB stars in Local Group and in Local Volume galax¬ 
ies (Bono. Salaris & Gilmozzil[201 3). 
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Fig. 1.— The Hertzsprung-Russell diagram for three sets of HB evolutionary models. From left to right black solid and dashed lines show the location of the 
ZAHB and of the central helium exhaustion, respectively. Red lines display HB evolutionary models populating the RRL instability strip. The individual stellar 
masses are 0.75, 0.76, 0.77 M/M Q for the most metal-poor (Z=0.0001, right panel), 0.65, 0.66, 0.67 M/M Q for the metal-intermediate (Z=0.001, middle panel) 
and 0.56, 0.57, 0.58 M/M Q for the more metal-rich (Z=0.008, left panel) chemical composition. 
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Fig. 2.— Location in the Hertzsprung-Russell diagram for a set of RRL models at fixed chemical composition (Z=0.0003, Y=0.245). The FU models are 
marked with filled circles, while the FOs with open circles. The black symbols mark pulsation models computed assuming the same stellar mass (0.716 M©) and 
three different luminosity levels: the Zero Age Horizontal Branch (ZAHB, sequence A), a luminosity level 0.1 dex brighter than the ZAHB (sequence B) and the 
luminosity level of central He exhaustion (sequence D). The red symbols display RRL models computed assuming a stellar mass ~10% smaller (0.65 M©) than 
the ZAHB mass value and 0.2 dex brighter than the ZAHB luminosity (sequence C). This sequence of pulsation models was computed to account for evolved 
RRLs. Similar sets of RRL models were computed for the other adopted chemical compositions (see Table ??). 
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Fig. 3.— Comparison between the FU periods (lo g(Pru)) predicted by individual models and the FU periods (log (Ppy )) given by the pulsation relation provided 
in this paper (black symbols). The red symbols display the same comparison, but using the pulsation relation provided by IDi Criscienzo, Marconi & Caputol 
( 223 ). 
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Fig. 4.— Relation between FU and FO periods for pulsation models located inside the so-called ”OR region”. The black solid line shows the linear regression 
based on theoretical models. The red and the green solid lines display the empirical and the classical relation between FU and FO periods (see text for more 
details). 
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Fig. 5.— Nonlinear total work integral as a function of integration time for three different FO (left panels) and three FU models (right panels). They are centrally 
located in the middle of the instability strip, and constructed assuming three different metal abundances, stellar masses and luminosity levels (see labeled values). 
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Fig. 6.— bolometric amplitude as a function of the integration time for the same models of Fig. \ 5 \ 
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Fig. 7.— From top to bottom the panels display predicted instability strips for FU (solid) and FO (dashed) pulsators at three different chemical compositions: 
Z = 0.0001 (top panel), Z = 0.001 (middle) and Z = 0.02 (bottom). The hot edges are plotted in blue, while the cool edges in red. The linear FOBE and FBE are 
represented by open and filled blue circles, respectively (see Appendix). 
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Fig. 8.— Comparison between the current instability strip boundaries (solid lines) and similar predictions bv IDi Criscienzo. Marconi & Caputol <20041) (dashed 
lines). The blue lines display the FO blue boundaries, while the red lines the FU red boundaries. 
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Fig. 9.— Tridimensional plot showing the correlation among stellar mass (log(M/M©)), metallicity (log(Z)) and period ratios (log (Pfo/Pfu)) f° r the pulsation 
models located inside the so-called ”OR region”. 
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Fig. 10.— Period-Radius relation for FU pulsators covering the entire range of chemical compositions adopted in the current investigation. The bl ack so lid line 
depicts the linear regression with a standard deviation of 0.03 dex. The red solid line shows the theoretical PR relation provided bv IMarco ni et all <20051), whil e 
the so lid green line the extrapolation to shorter periods of the PR for BL He rculis variables based o n pulsation models provided by Marconi & Di Criscienzo 
<20071) . The latter is almost identical to the empirical PR relation provided bv IBurki & Mevlar3 <19861) for the same class of variable stars and represented in the 
plot by the dashed green line. 
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Fig. 11.— Period-Radius-Metallicity (PRZ) relations for FU (top panel) and FO (bottom panel) pulsators with a standard deviation of 0.006 and 0.003 dex, 
respec tively. Black lines display the linear regression based on the current models, while the red lines show similar PRZ relations provided by IMarconi et all 
(2003). 
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Fig. 12.— Predicted bolometric light-curves for FU pulsators at fixed chemical composition (Z = 0.0006, Y = 0.245). The models refer to sequence (M 
0.67 Mq, logL/L© = 1.69). The effective temperatures and the pulsation periods are labeled. 
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Fig. 13. — Same as Figure IT2] but for FO pulsators. 
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Fig. 14.— Left panels — Predicted metal-dependent optical and NIR ( RIJHK ) PL relations FU and FO pulsators. Lines of different colors display predictions 
with metal abundances ranging from Z = 0.0001 (black lines) to Z = 0.02 (purple lines). Right panels - Same as the left ones, but for FU+FO pulsators. 
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Fig. 15.— Comparison among the current PLZ relations for the two different metal abundances (Z=0.0001, top; Z=0.001, bottom) and similar relations available 
in the literature. The standard deviations of the above PLZ relations, when available, are plotted in the top left comer of the top panel. 
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Fig. 16. — Left panels - Predicted metal-dependent optical PWZ relations. The color coding is the same as in Figure [14] The coefficients of the color terms are 
labeled on the Y-axis. Note either the minimal or the marginal metal dependence of the PW(V,B-V) (top panel) and of the PW(R,B-R) (second panel from top). 
Right panels - Same as the left ones, but for FU+FO pulsators. 
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Fig. 17.— Same as in Figure [I6l but for optical-NIR PWZ relations. 
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Fig. 18, 


Same as in Figure U6] but for NIR PWZ relations. 
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Fig. 19.— Comparison between the current PWZ(V,B-V) relations and similar relations provided bv IDi Criscienzo. Marconi & Canuto I <20041) for two different 
metal abundances: Z = 0.0001 (top panel) and Z = 0.001 (bottom panel). 
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Fig. 20. — Left Panels - Predicted metal-dependent triple bands optical-NIR PWZ relations. The color coding is the same as in Figure fl6l Right panels - Same 
as the left ones, but for FU+FO pulsators. 
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APPENDIX 
LINEAR MODELS 

The current set of linear models are purely radiative ( iBono & Stellingwerfil 1994 ) and do not provide an estimate of the red 
boundaries of the instability strip. The physical mechanism that is providing the quenching of radial oscillation is the increased 
efficiency of convective transport. However they provide the envelope structures and the linear radial eigenfunctions adopted by 
the nonlinear hydrodynamical models. The nonlinear analysis is pe rformed by imposing a constant perturbation velocity (5-10 
km se c -1 ) both to the FU and the FO linear radial eigenfunctions (Bon o & Stellingwerfil 19941 : IBono. Marconi. & S tellingwerf 
1999). After this initial perturbation the dynamical behavior is followed in time until the radial displacements approach their 
asymptotic behavior. In the current linear models the inner boundary of the static model is located at a distance of = 10% of the 
equilibrium photospheric radius. The outer boundary is located at an optical depth r = 0.0001 and the stellar mass attached to the 
outermost zone is at least one order of magnitude larger than the mass of the surface. The mass zoning of the stellar envelope was 
fixed assuming a mass ratio between consecutive zones of 1.04 for temperatures cooler than 60,000 K. It was increased by 0.001 
in the zones located at higher temperatures. As a whole we ended up with envelope models including typically 200-300 discrete 
mesh points (IBono et alJll998l) . allowing smooth variations of the physical parameters in both driving and damping envelope 
regions. The envelope mass decreases from 18% to a few percents of the total mass when moving from the blue faint to the red 
bright pulsation models. The dependence of the envelope mass on chemical composition is negligible. 

Fig. [7] shows the comparison between the nonlinear edges of the instability strip (§ 4) and the linear FU (filled circles) and 
FO (open circles) blue boundaries. We note that the linear FOBE and FBE shows a more evident dependence on the metal 
content than their nonlinear counterpart when moving from Z=0.0001 to Z-0.02. In particular, for the most metal-poor chemical 
composition (Z=0.0001, Y=0.245) the linear boundaries are significantly bluer than the nonlinear ones. This is an expected result, 
because the inclusion of convective transport causes a decrease in pulsation destabilization (driving), therefore, the blue edges 
move towards lower effective temperatures. The quenching caused by convection, close to the blue edges of the instability strip, 
becomes less efficient in the metal-r ich regime. It is partially counterbalanced by the driving of the K-bump, i.e. the radiative 
opacity peak located at =250,000 K (ISeaton et alJll994tlRogers & Iglesiaslll992tIBono et alJ[l996l) . 
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Table 1. Evolutionary lifetimes for typical ZAHB and MSTO stellar structures at different metal 

abundances. 


z a 

Y a 

T MS b 

W 

Thb/Tms 

M HB d 

Mro e 

0.0001 

0.245 

2056.59 

66.80 

0.03248 

0.80 

0.80 

0.0010 

0.246 

3433.67 

77.545 

0.02258 

0.65 

0.79 

0.0080 

0.256 

5969.38 

86.101 

0.01442 

0.57 

0.85 


a Metal (Z) and helium (Y) abundance by mass. 

b Central hydrogen burning lifetime between a point that is 0.25 dex 
fainter than the MSTO and the MSTO (Myr). 

c Central helium burning lifetime at 90% of core helium exhaustion 
(Myr). 

d ZAHB stellar mass roughly located at the center of the instability strip 
(solar units). 

e Stellar mass at the MSTO for an of 13 Gyrs (solar units). 
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Table 2. Intrinsic stellar parameters adopted to constract the pulsation models and predicted 
effective temperatures of the instability strip boundaries. The different sets of pulsation models 
were computed assuming the ZAHB mass value (column 1) and three different luminosity levels 
(column 2): the ZAHB level (sequence A), a luminosity level 0.1 dex brighter than the ZAHB 
(sequence B) and the luminosity level of central He exhaustion (sequence D). An additional 
sequence of RRL models was computed assuming a stellar mass ~10% smaller than the ZAHB 
mass value and 0.2 dex brighter than the ZAHB luminosity level (sequence C). 


Mzahb a 

log L/Lo b 

Sequence c 

FOBE i 

FBE i 

FORE a 

FRE a 



Z=0.0001 

Y=0.245 




0.80 

1.76 

A 

7250 

6850 

6550 

5950 

0.80 

1.86 

B 

7150 

6950 

6650 

5850 

0.80 

1.97 

D 

7050 

6850 

6650 

5750 

0.72 

1.96 

C 

6950 

6850 

6650 

5750 



Z=0.0003 

Y=0.245 




0.716 

1.72 

A 

7250 

6950 

6550 

5950 

0.716 

1.82 

B 

7150 

6950 

6650 

5850 

0.716 

1.99 

D 

6950 

6850 

6750 

5650 

0.65 

1.92 

C 

6950 

6850 

6350 

5750 



Z=0.0006 

Y=0.245 




0.67 

1.69 

A 

7250 

6950 

6650 

5950 

0.67 

1.79 

B 

7050 

6950 

6650 

5850 

0.67 

2.01 

D 


6850 


5750 

0.60 

1.89 

C 


6950 


5650 



Z=0.001 Y=0.246 




0.64 

1.67 

A 

7250 

7050 

6750 

5950 

0.64 

1.77 

B 

7050 

7050 

6750 

5850 

0.64 

1.99 

D 


6850 


5650 

0.58 

1.87 

C 


6950 


5850 



Z=0.004 Y=0.25 




0.59 

1.61 

A 

7250 

7050 

6750 

5950 

0.59 

1.71 

B 

7150 

6950 

6750 

5850 

0.59 

2.02 

D 


6750 


5650 

0.53 

1.81 

C 


6850 


5650 



Z=0.008 Y=0.256 




0.57 

1.58 

A 

7150 

7050 

6750 

5950 

0.57 

1.68 

B 

7050 

6950 

6150 

5850 

0.57 

2.02 

D 


6750 


5350 

0.51 

1.78 

C 


6850 


5550 



Z=0.02 Y=0.27 




0.54 

1.49 

A 

7150 

6950 

6750 

5950 

0.54 

1.59 

B 

7050 

6950 

6750 

5850 

0.54 

1.94 

D 


6650 


5250 

0.51 

1.69 

C 


6850 


5650 


a ZAHB stellar mass (solar units). 
b ogarithmic luminosity (solar units). 
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c Name of the sequence of pulsation models (see Figure ?? for more de¬ 
tails). 

d Effective temperature (K) of the instability strip boundaries: FOBE, FO 
blue edge; FBE, FU blue edge; FORE, FO red edge; FRE, FU red edge. 
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Table 3. Intensity-averaged mean magnitudes for the entire grid of FU models. For each set of 
models we list the chemical composition (metal and helium content), the stellar mass and the 
luminosity levels adopted for the different sequences. The first three columns give, for every 
pulsation model, the effective temperature, the logarithmic luminosity level and the logarithmic 
period. The last eight columns list the predicted optical Johnson-Cousins U, B, V, R, I and the NIR 
2MASS J, H, K intensity-averaged mean magnitudes. The last column gives the mean radius in 

solar units. 


T 

log L 

logP 

U 

B 

V 

R 

I 

J 

H 

K 

log (R/R 0 ) 

6800 

1.76 

-0.3016 

0.653 

Z= 

0.652 

=0.0001 Y: 

0.410 

=0.245 M: 

0.225 

=0.80 M 0 

-0.009 

-0.296 

-0.540 

-0.559 

0.7452 

6700 

1.76 

-0.2795 

0.672 

0.670 

0.407 

0.210 

-0.035 

-0.337 

-0.592 

-0.612 

0.7595 

6600 

1.76 

-0.2574 

0.691 

0.695 

0.409 

0.198 

-0.061 

-0.376 

-0.644 

-0.665 

0.7726 

6500 

1.76 

-0.2358 

0.710 

0.727 

0.416 

0.189 

-0.084 

-0.414 

-0.696 

-0.718 

0.7861 

6400 

1.76 

-0.2124 

0.731 

0.763 

0.426 

0.184 

-0.104 

-0.450 

-0.746 

-0.769 

0.7984 

6300 

1.76 

-0.1892 

0.750 

0.797 

0.438 

0.180 

-0.122 

-0.483 

-0.793 

-0.817 

0.8112 

6200 

1.76 

-0.1650 

0.768 

0.826 

0.447 

0.177 

-0.137 

-0.513 

-0.835 

-0.861 

0.8229 

6100 

1.76 

-0.1421 

0.786 

0.854 

0.457 

0.176 

-0.150 

-0.539 

-0.873 

-0.901 

0.8341 

6000 

1.76 

-0.1178 

0.802 

0.878 

0.465 

0.174 

-0.161 

-0.561 

-0.903 

-0.933 

0.8427 

6900 

1.86 

-0.2389 

0.385 

0.383 

0.164 

-0.007 

-0.228 

-0.502 

-0.735 

-0.754 

0.7816 

6800 

1.86 

-0.2180 

0.402 

0.397 

0.162 

-0.021 

-0.254 

-0.542 

-0.788 

-0.807 

0.7965 

6600 

1.86 

-0.1744 

0.440 

0.435 

0.158 

-0.050 

-0.308 

-0.624 

-0.893 

-0.914 

0.8228 

6400 

1.86 

-0.1282 

0.481 

0.500 

0.169 

-0.070 

-0.355 

-0.700 

-0.993 

-1.016 

0.8481 

6200 

1.86 

-0.0813 

0.521 

0.567 

0.191 

-0.077 

-0.390 

-0.764 

-1.084 

-1.110 

0.8728 

6100 

1.86 

-0.0567 

0.542 

0.599 

0.203 

-0.078 

-0.404 

-0.793 

-1.126 

-1.154 

0.8842 

6000 

1.86 

-0.0334 

0.563 

0.627 

0.213 

-0.079 

-0.416 

-0.819 

-1.164 

-1.194 

0.8958 

5900 

1.86 

-0.0088 

0.581 

0.651 

0.222 

-0.080 

-0.427 

-0.841 

-1.196 

-1.228 

0.9039 

6800 

1.97 

-0.1269 

0.123 

0.109 

-0.120 

-0.300 

-0.530 

-0.816 

-1.057 

-1.077 

0.8492 

6700 

1.97 

-0.1051 

0.143 

0.125 

-0.123 

-0.315 

-0.557 

-0.857 

-1.109 

-1.130 

0.8633 

6600 

1.97 

-0.0818 

0.164 

0.148 

-0.124 

-0.329 

-0.584 

-0.897 

-1.161 

-1.183 

0.8770 

6500 

1.97 

-0.0593 

0.186 

0.177 

-0.120 

-0.340 

-0.609 

-0.936 

-1.212 

-1.235 

0.8891 

6400 

1.97 

-0.0355 

0.207 

0.213 

-0.112 

-0.348 

-0.631 

-0.973 

-1.262 

-1.286 

0.9023 

6300 

1.97 

-0.0125 

0.228 

0.250 

-0.101 

-0.353 

-0.651 

-1.008 

-1.311 

-1.336 

0.9144 

6200 

1.97 

0.0131 

0.250 

0.285 

-0.089 

-0.356 

-0.668 

-1.042 

-1.358 

-1.385 

0.9279 

6100 

1.97 

0.0369 

0.273 

0.318 

-0.077 

-0.357 

-0.683 

-1.071 

-1.401 

-1.430 

0.9391 

6000 

1.97 

0.0623 

0.296 

0.348 

-0.066 

-0.358 

-0.696 

-1.100 

-1.443 

-1.474 

0.9519 

5900 

1.97 

0.0872 

0.318 

0.375 

-0.056 

-0.360 

-0.707 

-1.124 

-1.480 

-1.513 

0.9627 

5800 

1.97 

0.1114 

0.340 

0.401 

-0.046 

-0.360 

-0.718 

-1.146 

-1.514 

-1.549 

0.9713 

6800 

1.96 

-0.1034 

0.149 

Z: 

0.126 

=0.0001 Y: 

-0.103 

=0.245 M: 

-0.282 

=0.72 M 0 
-0.511 

-0.796 

-1.034 

-1.054 

0.8454 

6700 

1.96 

-0.0822 

0.168 

0.144 

-0.105 

-0.296 

-0.537 

-0.836 

-1.086 

-1.107 

0.8583 

6500 

1.96 

-0.0355 

0.212 

0.201 

-0.098 

-0.319 

-0.588 

-0.915 

-1.188 

-1.211 

0.8845 

6300 

1.96 

0.0112 

0.254 

0.273 

-0.078 

-0.330 

-0.628 

-0.985 

-1.286 

-1.311 

0.9095 

6100 

1.96 

0.0618 

0.299 

0.345 

-0.054 

-0.334 

-0.660 

-1.048 

-1.377 

-1.406 

0.9353 

6000 

1.96 

0.0865 

0.322 

0.379 

-0.043 

-0.335 

-0.672 

-1.075 

-1.418 

-1.449 

0.9469 

5800 

1.96 

0.1360 

0.369 

0.421 

-0.022 

-0.336 

-0.694 

-1.124 

-1.492 

-1.527 

0.9681 

6900 

1.72 

-0.3241 

0.743 

Z= 

0.732 

:0.0003 Y= 

0.504 

=0.245 M= 

0.330 

=0.716 M 0 

0.109 

-0.162 

-0.392 

-0.410 

0.7127 

6800 

1.72 

-0.3023 

0.761 

0.746 

0.500 

0.314 

0.082 

-0.203 

-0.445 

-0.464 

0.7260 

6700 

1.72 

-0.2811 

0.780 

0.765 

0.498 

0.300 

0.055 

-0.243 

-0.497 

-0.517 

0.7395 

6600 

1.72 

-0.2583 

0.801 

0.790 

0.499 

0.287 

0.029 

-0.283 

-0.549 

-0.570 

0.7531 

6500 

1.72 

-0.2363 

0.822 

0.822 

0.505 

0.277 

0.006 

-0.321 

-0.599 

-0.621 

0.7657 

6400 

1.72 

-0.2124 

0.842 

0.856 

0.513 

0.270 

-0.012 

-0.357 

-0.648 

-0.671 

0.7786 

6300 

1.72 

-0.1892 

0.861 

0.889 

0.523 

0.266 

-0.031 

-0.390 

-0.693 

-0.718 

0.7899 

6200 

1.72 

-0.1659 

0.882 

0.918 

0.533 

0.263 

-0.040 

-0.419 

-0.735 

-0.761 

0.8013 
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Table 3—Continued 


T f 

log L 

logP 

U 

B 

V 

R 

I 

J 

H 

K 

log (R/Ro) 

6100 

1.72 

-0.1420 

0.903 

0.948 

0.543 

0.261 

-0.064 

-0.446 

-0.775 

-0.803 

0.8135 

6000 

1.72 

-0.1175 

0.921 

0.973 

0.552 

0.260 

-0.079 

-0.469 

-0.806 

-0.837 

0.8227 

6900 

1.82 

-0.2414 

0.490 

0.475 

0.253 

0.081 

-0.138 

-0.410 

-0.638 

-0.657 

0.7621 

6800 

1.82 

-0.2192 

0.508 

0.489 

0.249 

0.065 

-0.166 

-0.450 

-0.691 

-0.711 

0.7752 

6600 

1.82 

-0.1747 

0.549 

0.529 

0.244 

0.035 

-0.227 

-0.532 

-0.794 

-0.816 

0.8018 

6400 

1.82 

-0.1287 

0.592 

0.595 

0.256 

0.016 

-0.267 

-0.607 

-0.893 

-0.917 

0.8275 

6200 

1.82 

-0.0806 

0.636 

0.663 

0.279 

0.009 

-0.300 

-0.671 

-0.985 

-1.012 

0.8520 

6000 

1.82 

-0.0326 

0.682 

0.722 

0.300 

0.006 

-0.322 

-0.726 

-1.065 

-1.096 

0.8740 

5900 

1.82 

-0.0088 

0.703 

0.747 

0.309 

0.005 

-0.330 

-0.749 

-1.098 

-1.131 

0.8837 

6800 

1.99 

-0.0771 

0.080 

0.040 

-0.192 

-0.371 

-0.597 

-0.877 

-1.110 

-1.130 

0.8585 

6700 

1.99 

-0.0539 

0.101 

0.061 

-0.192 

-0.384 

-0.623 

-0.917 

-1.163 

-1.184 

0.8728 

6600 

1.99 

-0.0313 

0.124 

0.086 

-0.191 

-0.397 

-0.649 

-0.958 

-1.215 

-1.237 

0.8861 

6500 

1.99 

-0.0077 

0.148 

0.120 

-0.186 

-0.407 

-0.673 

-0.996 

-1.266 

-1.289 

0.8986 

6400 

1.99 

0.0153 

0.172 

0.156 

-0.177 

-0.415 

-0.696 

-1.033 

-1.315 

-1.340 

0.9124 

6300 

1.99 

0.0399 

0.195 

0.192 

-0.166 

-0.419 

-0.715 

-1.068 

-1.364 

-1.390 

0.9248 

6200 

1.99 

0.0654 

0.219 

0.228 

-0.154 

-0.422 

-0.732 

-1.101 

-1.411 

-1.439 

0.9376 

6100 

1.99 

0.0895 

0.244 

0.260 

-0.143 

-0.424 

-0.746 

-1.131 

-1.456 

-1.485 

0.9504 

5900 

1.99 

0.1407 

0.297 

0.321 

-0.121 

-0.426 

-0.771 

-1.185 

-1.537 

-1.570 

0.9737 

5800 

1.99 

0.1645 

0.322 

0.348 

- 0.111 

-0.426 

-0.781 

-1.208 

-1.572 

-1.607 

0.9837 

5700 

1.99 

0.1903 

0.352 

0.377 

Z= 

-0.101 

=0.0003 Y= 

-0.426 

=0.245 M= 

-0.791 
=0.65 M 0 

-1.232 

-1.607 

-1.645 

0.9931 

6800 

1.92 

-0.107 

0.257 

0.222 

-0.015 

-0.195 

-0.423 

-0.704 

-0.938 

-0.958 

0.8239 

6700 

1.92 

-0.085 

0.278 

0.243 

-0.015 

-0.209 

-0.449 

-0.744 

-0.990 

-1.011 

0.8382 

6300 

1.92 

0.009 

0.367 

0.375 

0.013 

-0.242 

-0.539 

-0.892 

-1.189 

-1.215 

0.8899 

6100 

1.92 

0.060 

0.414 

0.439 

0.035 

-0.246 

-0.570 

-0.954 

-1.279 

-1.308 

0.9148 

5900 

1.92 

0.109 

0.466 

0.498 

0.057 

-0.248 

-0.593 

-1.007 

-1.358 

-1.391 

0.9365 

5800 

1.92 

0.134 

0.490 

0.525 

Z= 

0.066 

=0.0006 Y= 

-0.249 

=0.245 M= 

-0.604 

=0.67 M 0 

-1.030 

-1.392 

-1.427 

0.9463 

6900 

1.69 

-0.329 

0.825 

0.803 

0.571 

0.396 

0.177 

-0.092 

-0.321 

-0.339 

0.6971 

6800 

1.69 

-0.308 

0.843 

0.818 

0.566 

0.380 

0.150 

-0.133 

-0.374 

-0.393 

0.7115 

6700 

1.69 

-0.286 

0.865 

0.836 

0.564 

0.365 

0.123 

-0.174 

-0.426 

-0.446 

0.7245 

6600 

1.69 

-0.263 

0.887 

0.861 

0.565 

0.352 

0.097 

-0.214 

-0.477 

-0.498 

0.7370 

6500 

1.69 

-0.240 

0.909 

0.893 

0.569 

0.341 

0.072 

-0.252 

-0.527 

-0.549 

0.7507 

6400 

1.69 

-0.218 

0.930 

0.927 

0.576 

0.332 

0.050 

-0.289 

-0.575 

-0.599 

0.7628 

6300 

1.69 

-0.195 

0.952 

0.960 

0.587 

0.328 

0.031 

-0.322 

-0.621 

-0.646 

0.7752 

6200 

1.69 

-0.171 

0.974 

0.991 

0.596 

0.324 

0.016 

-0.352 

-0.664 

-0.691 

0.7863 

6100 

1.69 

-0.147 

0.998 

1.022 

0.607 

0.322 

0.002 

-0.380 

-0.704 

-0.733 

0.7987 

6000 

1.69 

-0.123 

1.020 

1.048 

0.615 

0.320 

-0.010 

-0.404 

-0.740 

-0.770 

0.8087 

6900 

1.79 

-0.245 

0.571 

0.541 

0.314 

0.142 

-0.074 

-0.341 

-0.566 

-0.585 

0.7470 

6800 

1.79 

-0.223 

0.590 

0.555 

0.310 

0.127 

-0.101 

-0.381 

-0.618 

-0.638 

0.7594 

6600 

1.79 

-0.179 

0.634 

0.598 

0.306 

0.097 

-0.155 

-0.463 

-0.721 

-0.743 

0.7868 

6400 

1.79 

-0.133 

0.680 

0.666 

0.320 

0.078 

-0.202 

-0.538 

-0.820 

-0.844 

0.8113 

6200 

1.79 

-0.084 

0.728 

0.734 

0.342 

0.071 

-0.237 

-0.603 

-0.913 

-0.940 

0.8365 

6000 

1.79 

-0.037 

0.777 

0.793 

0.361 

0.067 

-0.263 

-0.658 

-0.992 

-1.023 

0.8585 

5900 

1.79 

-0.012 

0.800 

0.820 

0.370 

0.065 

-0.275 

-0.682 

-1.026 

-1.060 

0.8694 

6800 

2.01 

-0.041 

0.050 

-0.012 

-0.280 

-0.470 

-0.691 

-0.958 

-1.164 

-1.184 

0.8647 
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Table 3—Continued 


Te 

log L 

logP 

U 

B 

V 

R 

I 

J 

H 

K 

log (R/R a ) 

6700 

2.01 

-0.017 

0.061 

0.000 

-0.258 

-0.450 

-0.685 

-0.976 

-1.216 

-1.238 

0.8822 

6600 

2.01 

0.006 

0.086 

0.028 

-0.255 

-0.461 

-0.710 

-1.016 

-1.269 

-1.292 

0.8959 

6400 

2.01 

0.054 

0.137 

0.099 

-0.241 

-0.479 

-0.757 

-1.092 

-1.370 

-1.396 

0.9232 

6300 

2.01 

0.078 

0.163 

0.135 

-0.230 

-0.484 

-0.776 

-1.127 

-1.419 

-1.446 

0.9355 

6200 

2.01 

0.103 

0.190 

0.172 

-0.218 

-0.487 

-0.794 

-1.160 

-1.467 

-1.495 

0.9483 

6000 

2.01 

0.154 

0.247 

0.238 

-0.196 

-0.491 

-0.823 

-1.221 

-1.556 

-1.588 

0.9737 

5900 

2.01 

0.180 

0.274 

0.268 

-0.186 

-0.493 

-0.834 

-1.247 

-1.593 

-1.627 

0.9844 

5800 

2.01 

0.210 

0.319 

0.309 

Z= 

-0.171 

=0.0006 Y: 

-0.493 

=0.245 M: 

-0.849 

=0.60 M 0 

-1.281 

-1.648 

-1.686 

0.986 

6900 

1.89 

-0.1319 

0.336 

0.272 

0.018 

-0.159 

-0.371 

-0.624 

-0.817 

-0.836 

0.7917 

6800 

1.89 

-0.1079 

0.342 

0.288 

0.047 

-0.134 

-0.359 

-0.637 

-0.867 

-0.887 

0.8096 

6700 

1.89 

-0.0859 

0.364 

0.310 

0.047 

-0.147 

-0.384 

-0.677 

-0.919 

-0.940 

0.8235 

6500 

1.89 

-0.0392 

0.413 

0.372 

0.055 

-0.170 

-0.435 

-0.755 

-1.020 

-1.044 

0.8492 

6300 

1.89 

0.0089 

0.461 

0.443 

0.075 

-0.181 

-0.474 

-0.825 

-1.118 

-1.144 

0.8753 

6000 

1.89 

0.0838 

0.541 

0.541 

0.108 

-0.188 

-0.519 

-0.916 

-1.251 

-1.282 

0.9121 

5900 

1.89 

0.1098 

0.568 

0.570 

0.117 

-0.189 

-0.531 

-0.942 

-1.290 

-1.323 

0.9220 

5800 

1.89 

0.1365 

0.599 

0.601 

0.128 

-0.190 

-0.542 

-0.967 

-1.327 

-1.363 

0.9340 

5700 

1.89 

0.1610 

0.642 

0.637 

Z: 

0.141 

=0.001 Y= 

-0.190 

=0.245 M= 

-0.554 

=0.64 M 0 

-0.997 

-1.375 

-1.413 

0.9468 

7000 

1.67 

-0.3518 

0.865 

0.837 

0.617 

0.451 

0.245 

-0.009 

-0.223 

-0.240 

0.6736 

6900 

1.67 

-0.3304 

0.881 

0.850 

0.613 

0.436 

0.219 

-0.049 

-0.275 

-0.293 

0.6875 

6800 

1.67 

-0.3090 

0.902 

0.866 

0.608 

0.420 

0.191 

-0.090 

-0.327 

-0.346 

0.7017 

6700 

1.67 

-0.2872 

0.923 

0.885 

0.605 

0.405 

0.164 

-0.130 

-0.378 

-0.398 

0.7137 

6600 

1.67 

-0.2651 

0.945 

0.909 

0.605 

0.391 

0.137 

-0.170 

-0.429 

-0.450 

0.7265 

6500 

1.67 

-0.2421 

0.967 

0.942 

0.610 

0.380 

0.113 

-0.208 

-0.477 

-0.500 

0.7391 

6400 

1.67 

-0.2191 

0.990 

0.976 

0.619 

0.373 

0.092 

-0.243 

-0.525 

-0.549 

0.7517 

6300 

1.67 

-0.1954 

1.013 

1.009 

0.629 

0.369 

0.075 

-0.276 

-0.571 

-0.596 

0.7639 

6200 

1.67 

-0.1727 

1.037 

1.040 

0.639 

0.366 

0.059 

-0.305 

-0.613 

-0.640 

0.7765 

6100 

1.67 

-0.1482 

1.061 

1.069 

0.648 

0.364 

0.046 

-0.332 

-0.652 

-0.681 

0.7876 

6000 

1.67 

-0.1247 

1.083 

1.095 

0.656 

0.361 

0.034 

-0.356 

-0.687 

-0.717 

0.7972 

7000 

1.77 

-0.2614 

0.616 

0.574 

0.355 

0.190 

-0.016 

-0.271 

-0.483 

-0.501 

0.7272 

6900 

1.77 

-0.2395 

0.634 

0.590 

0.353 

0.177 

-0.040 

-0.309 

-0.533 

-0.552 

0.7406 

6800 

1.77 

-0.2178 

0.655 

0.606 

0.349 

0.161 

-0.067 

-0.349 

-0.584 

-0.604 

0.7537 

6700 

1.77 

-0.1974 

0.676 

0.626 

0.346 

0.146 

-0.094 

-0.389 

-0.634 

-0.655 

0.7661 

6500 

1.77 

-0.1515 

0.723 

0.690 

0.356 

0.125 

-0.142 

-0.465 

-0.733 

-0.756 

0.7912 

6300 

1.77 

-0.1052 

0.771 

0.759 

0.376 

0.116 

-0.180 

-0.531 

-0.826 

-0.852 

0.8154 

6100 

1.77 

-0.0573 

0.820 

0.818 

0.396 

0.111 

-0.208 

-0.588 

-0.909 

-0.938 

0.8384 

6000 

1.77 

-0.0340 

0.845 

0.846 

0.405 

0.109 

-0.220 

-0.613 

-0.945 

-0.976 

0.8499 

5900 

1.77 

-0.0100 

0.869 

0.872 

0.413 

0.107 

-0.231 

-0.636 

-0.978 

-1.011 

0.8582 

6800 

1.99 

-0.0410 

0.097 

0.023 

-0.219 

-0.398 

-0.618 

-0.892 

-1.114 

-1.135 

0.8582 

6700 

1.99 

-0.0185 

0.122 

0.049 

-0.216 

-0.408 

-0.642 

-0.932 

-1.168 

-1.190 

0.8731 

6500 

1.99 

0.0290 

0.174 

0.114 

-0.206 

-0.429 

-0.691 

-1.010 

-1.271 

-1.295 

0.8992 

6400 

1.99 

0.0526 

0.200 

0.151 

-0.196 

-0.436 

-0.712 

-1.046 

-1.321 

-1.346 

0.9121 

6300 

1.99 

0.0770 

0.227 

0.188 

-0.185 

-0.441 

-0.731 

-1.081 

-1.370 

-1.397 

0.9256 

6000 

1.99 

0.1547 

0.310 

0.285 

-0.153 

-0.448 

-0.775 

-1.169 

-1.498 

-1.530 

0.9611 

5900 

1.99 

0.1839 

0.347 

0.321 

-0.141 

-0.449 

-0.789 

-1.200 

-1.545 

-1.579 

0.9745 
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Table 3—Continued 


T f 

log L 

logf 

U 

B 

V 

R 

I 

J 

H 

K 

log (R/Ro) 

5800 

1.99 

0.2091 

0.395 

0.362 

-0.127 

-0.449 

-0.803 

-1.233 

-1.599 

-1.635 

0.9899 

5700 

1.99 

0.2320 

0.436 

0.394 

Z= 

-0.116 

=0.001 Y= 

-0.448 

=0.245 M= 

-0.813 

=0.58 Mq 

-1.259 

-1.641 

-1.679 

1.0021 

6900 

1.87 

-0.1259 

0.386 

0.319 

0.085 

-0.088 

-0.302 

-0.568 

-0.603 

-0.803 

0.7903 

6700 

1.87 

-0.0819 

0.431 

0.367 

0.088 

-0.111 

-0.351 

-0.645 

-0.886 

-0.908 

0.8165 

6500 

1.87 

-0.0355 

0.480 

0.434 

0.101 

-0.130 

-0.397 

-0.720 

-0.985 

-1.010 

0.8422 

6400 

1.87 

-0.0126 

0.504 

0.468 

0.111 

-0.135 

-0.416 

-0.754 

-1.033 

-1.058 

0.8552 

6300 

1.87 

0.0120 

0.529 

0.502 

0.121 

-0.139 

-0.434 

-0.786 

-1.079 

-1.105 

0.8664 

6100 

1.87 

0.0606 

0.582 

0.565 

0.142 

-0.144 

-0.463 

-0.845 

-1.165 

-1.195 

0.8904 

5900 

1.87 

0.1109 

0.637 

0.623 

Z: 

0.161 

=0.004 Y= 

-0.147 

=0.25 M= 

-0.486 

=0.59 M 0 

-0.895 

-1.240 

-1.272 

0.9115 

7000 

1.61 

-0.3608 

1.041 

0.971 

0.728 

0.558 

0.356 

0.111 

-0.092 

0.110 

0.6463 

6900 

1.61 

-0.3404 

1.062 

0.987 

0.725 

0.544 

0.330 

0.072 

-0.142 

-0.161 

0.6595 

6800 

1.61 

-0.3215 

1.084 

1.004 

0.722 

0.530 

0.307 

0.037 

-0.186 

-0.206 

0.6700 

6700 

1.61 

-0.2974 

1.113 

1.032 

0.722 

0.515 

0.278 

-0.006 

-0.241 

-0.263 

0.6845 

6600 

1.61 

-0.2800 

1.134 

1.056 

0.723 

0.504 

0.256 

-0.039 

-0.283 

-0.306 

0.6952 

6500 

1.61 

-0.2539 

1.166 

1.097 

0.731 

0.494 

0.230 

-0.081 

-0.340 

-0.364 

0.7103 

6400 

1.61 

-0.2312 

1.194 

1.131 

0.740 

0.487 

0.210 

-0.116 

-0.387 

-0.412 

0.7228 

6300 

1.61 

-0.2088 

1.223 

1.165 

0.750 

0.483 

0.192 

-0.148 

-0.433 

-0.459 

0.7347 

6200 

1.61 

-0.1862 

1.252 

1.195 

0.758 

0.479 

0.177 

-0.177 

-0.473 

-0.501 

0.7461 

6100 

1.61 

-0.1648 

1.282 

1.224 

0.766 

0.476 

0.164 

-0.203 

-0.511 

-0.541 

0.7569 

6000 

1.61 

-0.1412 

1.314 

1.254 

0.774 

0.474 

0.151 

-0.230 

-0.549 

-0.581 

0.7682 

6900 

1.71 

-0.2576 

0.812 

0.727 

0.469 

0.291 

0.079 

-0.178 

-0.388 

-0.408 

0.7084 

6800 

1.71 

-0.2353 

0.838 

0.749 

0.469 

0.277 

0.054 

-0.217 

-0.439 

-0.460 

0.7210 

6700 

1.71 

-0.2193 

0.859 

0.768 

0.468 

0.266 

0.033 

-0.249 

-0.480 

-0.502 

0.7324 

6500 

1.71 

-0.1694 

0.921 

0.841 

0.478 

0.242 

-0.020 

-0.331 

-0.588 

-0.612 

0.7598 

6300 

1.71 

-0.1240 

0.980 

0.911 

0.497 

0.231 

-0.058 

-0.399 

-0.683 

-0.710 

0.7847 

6100 

1.71 

-0.0764 

1.046 

0.977 

0.516 

0.225 

-0.089 

-0.459 

-0.770 

-0.800 

0.8087 

5900 

1.71 

-0.0320 

1.114 

1.039 

0.534 

0.221 

-0.114 

-0.510 

-0.844 

-0.878 

0.8300 

6700 

2.02 

0.0481 

0.096 

-0.026 

-0.323 

-0.522 

-0.752 

-1.034 

-1.258 

-1.281 

0.8887 

6300 

2.02 

0.1509 

0.231 

0.123 

-0.287 

-0.550 

-0.838 

-1.184 

-1.467 

-1.496 

0.9452 

6200 

2.02 

0.1648 

0.252 

0.144 

-0.281 

-0.553 

-0.849 

-1.204 

-1.496 

-1.526 

0.9519 

6100 

2.02 

0.2056 

0.323 

0.210 

-0.261 

-0.557 

-0.875 

-1.258 

-1.578 

-1.610 

0.9755 

5900 

2.02 

0.2547 

0.433 

0.297 

-0.232 

-0.558 

-0.904 

-1.326 

-1.684 

-1.719 

1.0047 

5800 

2.02 

0.2792 

0.482 

0.334 

-0.220 

-0.557 

-0.914 

-1.351 

-1.723 

-1.761 

1.0159 

5700 

2.02 

0.3017 

0.524 

0.366 

Z: 

-0.210 

=0.004 Y: 

-0.557 

=0.25 M= 

-0.924 

=0.53 M 0 

-1.374 

-1.756 

-1.796 

1.0266 

6800 

1.81 

-0.1207 

0.590 

0.483 

0.207 

0.019 

-0.202 

-0.471 

-0.689 

-0.711 

0.7719 

6700 

1.81 

-0.1015 

0.616 

0.509 

0.210 

0.009 

-0.223 

-0.505 

-0.734 

-0.757 

0.7832 

6500 

1.81 

-0.0535 

0.678 

0.582 

0.223 

-0.011 

-0.272 

-0.584 

-0.840 

-0.865 

0.8118 

6400 

1.81 

-0.0303 

0.709 

0.619 

0.232 

-0.017 

-0.292 

-0.619 

-0.888 

-0.914 

0.8234 

6100 

1.81 

0.0387 

0.804 

0.718 

0.260 

-0.029 

-0.340 

-0.710 

-1.018 

-1.049 

0.8599 

6000 

1.81 

0.0621 

0.839 

0.752 

0.270 

-0.031 

-0.354 

-0.738 

-1.058 

-1.091 

0.8712 

5800 

1.81 

0.1195 

0.951 

0.840 

0.298 

-0.033 

-0.384 

-0.807 

-1.164 

-1.201 

0.9011 

5700 

1.81 

0.1449 

1.000 

0.874 

Z= 

0.309 

=0.008 Y= 

-0.032 

=0.256 M= 

-0.393 

=0.57 M 0 

-0.832 

-1.204 

-1.243 

0.9132 
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Table 3—Continued 


T f 

log L 

logP 

U 

B 

V 

R 

I 

J 

H 

K 

log (R/Ro) 

7000 

1.58 

-0.3677 

1.128 

1.026 

0.767 

0.595 

0.397 

0.159 

-0.033 

-0.052 

0.6326 

6900 

1.58 

-0.3465 

1.152 

1.047 

0.767 

0.583 

0.373 

0.122 

-0.083 

-0.103 

0.6459 

6800 

1.58 

-0.3251 

1.179 

1.068 

0.765 

0.569 

0.347 

0.083 

-0.132 

-0.153 

0.6583 

6600 

1.58 

-0.2826 

1.240 

1.128 

0.769 

0.544 

0.297 

0.007 

-0.230 

-0.254 

0.6839 

6400 

1.58 

-0.2365 

1.304 

1.200 

0.784 

0.528 

0.253 

-0.065 

-0.328 

-0.354 

0.7089 

6200 

1.58 

-0.1900 

1.375 

1.268 

0.802 

0.519 

0.220 

-0.128 

-0.416 

-0.445 

0.7344 

6000 

1.58 

-0.1458 

1.455 

1.337 

0.821 

0.513 

0.191 

-0.185 

-0.499 

-0.532 

0.7565 

6900 

1.68 

-0.2668 

0.914 

0.798 

0.522 

0.340 

0.133 

-0.117 

-0.319 

-0.339 

0.6934 

6800 

1.68 

-0.2462 

0.942 

0.821 

0.522 

0.328 

0.108 

-0.155 

-0.368 

-0.390 

0.7066 

6600 

1.68 

-0.2024 

1.004 

0.883 

0.527 

0.304 

0.058 

-0.232 

-0.468 

-0.492 

0.7324 

6400 

1.68 

-0.1574 

1.070 

0.954 

0.542 

0.288 

0.015 

-0.303 

-0.564 

-0.591 

0.7579 

6200 

1.68 

-0.1115 

1.144 

1.026 

0.561 

0.279 

-0.020 

-0.368 

-0.657 

-0.686 

0.7818 

6000 

1.68 

-0.0642 

1.228 

1.099 

0.581 

0.273 

-0.050 

-0.428 

-0.743 

-0.776 

0.8059 

5900 

1.68 

-0.0402 

1.276 

1.136 

0.592 

0.271 

-0.063 

-0.456 

-0.784 

-0.819 

0.8177 

6680 

2.02 

0.0683 

0.136 

-0.008 

-0.346 

-0.557 

-0.789 

-1.061 

-1.270 

-1.295 

0.8892 

6300 

2.02 

0.1645 

0.283 

0.120 

-0.306 

-0.569 

-0.853 

1.193 

-1.473 

-1.503 

0.9466 

6200 

2.02 

0.1911 

0.334 

0.163 

-0.294 

-0.573 

-0.872 

-1.231 

-1.529 

-1.560 

0.9606 

5800 

2.02 

0.2892 

0.588 

0.362 

-0.233 

-0.578 

-0.937 

-1.374 

-1.742 

-1.781 

1.0203 

5600 

2.02 

0.3390 

0.717 

0.453 

-0.205 

-0.577 

-0.959 

-1.428 

-1.824 

-1.867 

1.0450 

5400 

2.02 

0.3871 

0.850 

0.538 

Z= 

-0.178 

0.008 Y: 

-0.573 

=0.256 M 

-0.978 

=0.51 M 0 

-1.476 

1.990 

-1.943 

1.0656 

6800 

1.78 

-0.1299 

0.695 

0.558 

0.263 

0.072 

-0.146 

-0.408 

-0.617 

-0.640 

0.7561 

6700 

1.78 

-0.1082 

0.727 

0.588 

0.268 

0.061 

-0.169 

-0.446 

-0.668 

-0.692 

0.7703 

6500 

1.78 

-0.0621 

0.795 

0.660 

0.280 

0.042 

-0.215 

-0.521 

-0.768 

-0.795 

0.7956 

6100 

1.78 

0.0314 

0.948 

0.809 

0.318 

0.022 

-0.288 

-0.653 

-0.953 

-0.986 

0.8466 

5900 

1.78 

0.0851 

1.050 

0.890 

0.341 

0.018 

-0.318 

-0.715 

-1.045 

-1.082 

0.8724 

5600 

1.78 

0.1594 

1.259 

1.029 

Z 

0.388 

=0.02 Y: 

0.023 

=0.27 M= 

-0.353 

0.54 M g 

-0.810 

-1.195 

-1.238 

0.9156 

6900 

1.49 

-0.4014 

1.439 

1.283 

0.968 

0.777 

0.576 

0.336 

0.146 

0.125 

0.5960 

6800 

1.49 

-0.3795 

1.473 

1.311 

0.969 

0.766 

0.553 

0.298 

0.098 

0.074 

0.6091 

6700 

1.49 

-0.3557 

1.511 

1.344 

0.972 

0.752 

0.526 

0.258 

0.045 

0.021 

0.6228 

6500 

1.49 

-0.3129 

1.587 

1.416 

0.984 

0.734 

0.482 

0.186 

-0.050 

-0.077 

0.6494 

6200 

1.49 

-0.2455 

1.721 

1.528 

1.010 

0.718 

0.426 

0.087 

-0.188 

-0.219 

0.6871 

6000 

1.49 

-0.2008 

1.829 

1.607 

1.031 

0.712 

0.396 

0.026 

-0.274 

-0.309 

0.7105 

6900 

1.59 

-0.3148 

1.194 

1.026 

0.711 

0.521 

0.322 

0.082 

-0.104 

-0.126 

0.6471 

6700 

1.59 

-0.2725 

1.264 

1.086 

0.719 

0.502 

0.277 

0.009 

-0.203 

-0.227 

0.6735 

6500 

1.59 

-0.2295 

1.340 

1.158 

0.731 

0.484 

0.233 

-0.063 

-0.298 

-0.325 

0.6988 

6300 

1.59 

-0.1834 

1.430 

1.236 

0.749 

0.471 

0.193 

-0.133 

-0.394 

-0.424 

0.7242 

6100 

1.59 

-0.1372 

1.534 

1.317 

0.769 

0.463 

0.159 

-0.197 

-0.485 

-0.518 

0.7495 

6000 

1.59 

-0.1136 

1.593 

1.358 

0.780 

0.461 

0.144 

-0.228 

-0.529 

-0.564 

0.7632 

5900 

1.59 

-0.0943 

1.647 

1.395 

0.790 

0.459 

0.132 

-0.254 

-0.566 

-0.603 

0.7727 

6600 

1.94 

0.0474 

0.440 

0.238 

-0.164 

-0.397 

-0.632 

-0.911 

-1.121 

-1.148 

0.8604 

6400 

1.94 

0.0940 

0.530 

0.303 

-0.144 

-0.401 

-0.662 

-0.975 

-1.222 

-1.252 

0.8894 

6200 

1.94 

0.1410 

0.633 

0.379 

-0.120 

-0.405 

-0.692 

-1.042 

-1.323 

-1.356 

0.9190 

6000 

1.94 

0.1960 

0.786 

0.489 

-0.087 

-0.408 

-0.728 

-1.124 

-1.445 

-1.482 

0.9526 

5700 

1.94 

0.2685 

1.025 

0.652 

-0.034 

-0.406 

-0.768 

-1.226 

-1.596 

-1.639 

0.9967 
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Table 3—Continued 


T e 

log L 

logP 

U 

B 

V 

R 

I 

J 

H 

K 

log (R/R q ) 

5400 

1.94 

0.3444 

1.297 

0.814 

0.019 

-0.398 

-0.797 

-1.310 

-1.720 

-1.768 

1.0321 

5300 

1.94 

0.3699 

1.397 

0.870 

0.039 

-0.393 

-0.803 

-1.335 

-1.758 

-1.809 

1.0449 

5200 

1.94 

0.3962 

1.505 

0.923 

Z: 

0.053 

=0.02 Y= 

-0.392 

=0.27 M= 

-0.813 

=0.51 M g 

-1.357 

-1.786 

-1.839 

1.0518 

6800 

1.69 

-0.1933 

0.981 

0.799 

0.457 

0.255 

0.045 

-0.207 

-0.400 

-0.424 

0.7098 

6700 

1.69 

-0.1715 

1.016 

0.828 

0.465 

0.250 

0.027 

-0.241 

-0.451 

-0.476 

0.7231 

6500 

1.69 

-0.1268 

1.098 

0.902 

0.478 

0.233 

-0.017 

-0.314 

-0.550 

-0.577 

0.7502 

6200 

1.69 

-0.0580 

1.242 

1.024 

0.507 

0.215 

-0.075 

-0.418 

-0.693 

-0.725 

0.7886 

5900 

1.69 

0.0122 

1.423 

1.153 

0.541 

0.207 

-0.122 

-0.512 

-0.826 

-0.864 

0.8270 

5700 

1.69 

0.0638 

1.580 

1.251 

0.572 

0.208 

-0.147 

-0.576 

-0.922 

-0.963 

0.8549 







Te 

7200 

7100 

7000 

6900 

6800 

6700 

6600 

7100 

7000 

6900 

6800 

6700 

7000 

6900 

6800 

6700 

6900 

6700 

7200 

7100 

7000 

6900 

6800 

6700 

6600 

7100 

7000 

6900 

6800 

6700 

6900 

6800 

6700 

6900 

6700 

6600 

6500 

6400 

7200 

7100 

7000 
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Table 4. Same as Tabled but for FO models. 


log L 

logP 

U 

B 

V 

R 

I 

J 

H 

K 

log(R/R e ) 

1.76 

-0.5088 

0.613 

Z: 

0.562 

=0.0001 Y= 

0.366 

=0.245 M= 

0.223 

=0.72 M g 

0.040 

-0.182 

-0.351 

-0.367 

0.6907 

1.76 

-0.4901 

0.620 

0.583 

0.378 

0.226 

0.030 

-0.209 

-0.398 

-0.443 

0.7039 

1.76 

-0.4703 

0.632 

0.604 

0.382 

0.218 

0.010 

-0.243 

-0.446 

-0.463 

0.7178 

1.76 

-0.4495 

0.645 

0.625 

0.385 

0.209 

-0.012 

-0.279 

-0.494 

-0.512 

0.7299 

1.76 

-0.4282 

0.660 

0.652 

0.389 

0.198 

-0.036 

-0.316 

-0.544 

-0.563 

0.726 

1.76 

-0.4077 

0.675 

0.682 

0.396 

0.191 

-0.058 

-0.352 

-0.592 

-0.612 

0.7548 

1.76 

-0.3855 

0.689 

0.711 

0.405 

0.186 

-0.076 

-0.385 

-0.638 

-0.659 

0.7670 

1.86 

-0.4093 

0.375 

0.325 

0.109 

-0.045 

-0.239 

-0.472 

-0.646 

-0.663 

0.7526 

1.86 

-0.3891 

0.384 

0.343 

0.117 

-0.046 

-0.252 

-0.501 

-0.693 

-0.711 

0.7651 

1.86 

-0.3691 

0.397 

0.367 

0.124 

-0.052 

-0.271 

-0.534 

-0.741 

-0.760 

0.7794 

1.86 

-0.3472 

0.411 

0.394 

0.131 

-0.058 

-0.290 

-0.568 

-0.791 

-0.810 

0.7926 

1.86 

-0.3261 

0.426 

0.421 

0.140 

-0.063 

-0.309 

-0.602 

-0.839 

-0.859 

0.8057 

1.97 

-0.3011 

0.111 

0.063 

-0.168 

-0.332 

-0.538 

-0.782 

-0.966 

-0.983 

0.8190 

1.97 

-0.2797 

0.123 

0.082 

-0.159 

-0.333 

-0.550 

-0.810 

-1.013 

-1.031 

0.8330 

1.97 

-0.2583 

0.137 

0.107 

-0.151 

-0.337 

-0.567 

-0.843 

-1.062 

-1.081 

0.8466 

1.97 

-0.2368 

0.152 

0.134 

-0.142 

-0.342 

-0.586 

-0.877 

-1.110 

-1.131 

0.8596 

1.96 

-0.2596 

0.148 

Z: 

0.105 

=0.0001 Y= 

-0.137 

=0.245 M= 

-0.311 

=0.72 M 0 

-0.528 

-0.788 

-1.001 

-1.006 

0.8279 

1.96 

-0.2157 

0.176 

0.154 

-0.120 

-0.319 

-0.562 

-0.853 

-1.098 

-1.108 

0.8545 

1.72 

-0.5123 

0.718 

Z= 

0.657 

=0.0003 Y= 

0.452 

:0.245 M= 

0.309 

:0.716 M 0 

0.127 

-0.091 

-0.254 

-0.2705 

0.6704 

1.72 

-0.4926 

0.726 

0.675 

0.462 

0.308 

0.114 

-0.119 

-0.302 

-0.318 

0.6830 

1.72 

-0.4732 

0.739 

0.696 

0.466 

0.301 

0.094 

-0.153 

-0.349 

-0.366 

0.6960 

1.72 

-0.4516 

0.753 

0.721 

0.471 

0.293 

0.073 

-0.188 

-0.397 

-0.415 

0.7089 

1.72 

-0.4305 

0.768 

0.748 

0.477 

0.285 

0.052 

-0.224 

-0.446 

-0.465 

0.7230 

1.72 

-0.4098 

0.783 

0.777 

0.485 

0.278 

0.032 

-0.258 

-0.494 

-0.514 

0.7350 

1.72 

-0.3881 

0.801 

0.808 

0.496 

0.276 

0.015 

-0.289 

-0.538 

-0.560 

0.7468 

1.82 

-0.4121 

0.478 

0.418 

0.198 

0.044 

-0.149 

-0.377 

-0.549 

-0.566 

0.7321 

1.82 

-0.3921 

0.489 

0.438 

0.206 

0.041 

-0.163 

-0.407 

-0.595 

-0.613 

0.7445 

1.82 

-0.3709 

0.504 

0.462 

0.213 

0.036 

-0.181 

-0.440 

-0.643 

-0.662 

0.7578 

1.82 

-0.3498 

0.518 

0.487 

0.220 

0.030 

-0.200 

-0.474 

-0.692 

-0.712 

0.7710 

1.82 

-0.3295 

0.534 

0.515 

0.229 

0.025 

-0.219 

-0.508 

-0.740 

-0.761 

0.7842 

1.99 

-0.2342 

0.078 

0.024 

-0.225 

-0.400 

-0.615 

-0.870 

-1.064 

-1.083 

0.8406 

1.99 

-0.2104 

0.092 

0.044 

-0.215 

-0.401 

-0.628 

-0.900 

-1.113 

-1.133 

0.8564 

1.99 

-0.1891 

0.109 

0.069 

-0.206 

-0.404 

-0.645 

-0.933 

-1.164 

-1.185 

0.8706 

1.92 

-0.2644 

0.253 

Z: 

0.205 

=0.0003 Y= 

-0.045 

=0.245 M= 

-0.221 

=0.65 M 0 

-0.437 

-0.693 

-0.890 

-0.909 

0.8073 

1.92 

-0.2194 

0.284 

0.252 

-0.026 

-0.226 

-0.468 

-0.758 

-0.990 

-1.011 

0.8359 

1.92 

-0.1956 

0.303 

0.283 

-0.017 

-0.231 

-0.488 

-0.794 

-1.043 

-1.065 

0.8501 

1.92 

-0.1731 

0.324 

0.316 

-0.007 

-0.236 

-0.508 

-0.829 

-1.095 

-1.118 

0.8638 

1.92 

-0.1515 

0.345 

0.349 

0.004 

-0.240 

-0.526 

-0.864 

-1.145 

-1.169 

0.8762 

1.69 

-0.5192 

0.798 

Z: 

0.730 

=0.0006 Y= 

0.515 

=0.245 M= 

0.370 

=0.67 M 0 
0.190 

-0.025 

-0.183 

-0.199 

0.6548 

1.69 

-0.4991 

0.807 

0.746 

0.524 

0.369 

0.178 

-0.053 

-0.230 

-0.246 

0.6684 

1.69 

-0.4787 

0.820 

0.768 

0.530 

0.363 

0.159 

-0.086 

-0.278 

-0.295 

0.6815 





Table 4—Continued 


T ( , 

log L 

logP 

U 

B 

V 

R 

I 

J 

H 

K 

log(R/R e ) 

6900 

1.69 

-0.4569 

0.835 

0.792 

0.535 

0.355 

0.138 

-0.121 

-0.326 

-0.344 

0.6945 

6800 

1.69 

-0.4362 

0.850 

0.819 

0.541 

0.348 

0.118 

-0.155 

-0.374 

-0.393 

0.7076 

6700 

1.69 

-0.4162 

0.867 

0.847 

0.549 

0.342 

0.098 

-0.189 

-0.422 

-0.442 

0.7193 

7000 

1.79 

-0.3982 

0.570 

0.509 

0.271 

0.105 

-0.097 

-0.339 

-0.524 

-0.542 

0.7293 

6900 

1.79 

-0.3757 

0.584 

0.532 

0.277 

0.100 

-0.115 

-0.371 

-0.572 

-0.591 

0.7431 

6800 

1.79 

-0.3547 

0.600 

0.558 

0.284 

0.094 

-0.134 

-0.405 

-0.621 

-0.641 

0.7569 

6700 

1.79 

-0.3338 

0.618 

0.585 

Z= 

0.292 

0.001 Y: 

0.088 

=0.245 M 

-0.153 
=0.64 M 0 

-0.439 

-0.670 

-0.691 

0.7703 

7200 

1.67 

-0.5218 

0.850 

0.774 

0.560 

0.415 

0.237 

0.023 

-0.135 

-0.151 

0.6441 

7000 

1.67 

-0.4813 

0.875 

0.816 

0.573 

0.405 

0.202 

-0.041 

-0.229 

-0.247 

0.6704 

6900 

1.67 

-0.4602 

0.890 

0.840 

0.578 

0.398 

0.182 

-0.075 

-0.278 

-0.296 

0.6845 

6800 

1.67 

-0.4397 

0.907 

0.867 

0.584 

0.390 

0.161 

-0.109 

-0.326 

-0.345 

0.6964 

7000 

1.77 

-0.3931 

0.629 

0.563 

0.317 

0.148 

-0.057 

-0.302 

-0.490 

-0.508 

0.7232 

6900 

1.77 

-0.3720 

0.645 

0.587 

0.323 

0.141 

-0.075 

-0.334 

-0.537 

-0.556 

0.7368 

6800 

1.77 

-0.3505 

0.662 

0.613 

Z= 

0.330 

0.004 Y= 

0.136 

=0.245 M: 

-0.094 

=0.59 M 0 

-0.368 

-0.585 

-0.605 

0.7504 

7200 

1.61 

-0.5398 

1.021 

0.922 

0.684 

0.533 

0.357 

0.149 

-0.004 

-0.020 

0.6162 

6900 

1.61 

-0.4775 

1.070 

0.988 

0.701 

0.515 

0.302 

0.052 

-0.144 

-0.163 

0.6563 

6800 

1.61 

-0.4601 

1.088 

1.011 

0.706 

0.508 

0.284 

0.022 

-0.186 

-0.206 

0.6676 

7100 

1.71 

-0.4361 

0.789 

0.688 

0.434 

0.272 

0.085 

-0.135 

-0.298 

-0.316 

0.6788 

7000 

1.71 

-0.4161 

0.804 

0.706 

0.441 

0.269 

0.071 

-0.164 

-0.345 

-0.363 

0.6930 

6900 

1.71 

-0.3953 

0.822 

0.729 

0.447 

0.263 

0.052 

-0.197 

-0.392 

-0.412 

0.7052 

6800 

1.71 

-0.3739 

0.843 

0.756 

Z= 

0.454 

0.008 Y= 

0.256 

=0.245 M: 

0.033 

=0.57 M 0 

-0.231 

-0.440 

-0.461 

0.7195 

7100 

1.58 

-0.5303 

1.118 

1.001 

0.737 

0.572 

0.388 

0.172 

0.009 

-0.009 

0.6174 

7000 

1.58 

-0.5085 

1.138 

1.025 

0.742 

0.565 

0.369 

0.139 

-0.037 

-0.056 

0.6303 

6900 

1.58 

-0.4877 

1.159 

1.050 

0.747 

0.557 

0.348 

0.105 

-0.084 

-0.104 

0.6420 

6800 

1.58 

-0.4666 

1.182 

1.078 

0.752 

0.549 

0.328 

0.071 

-0.131 

-0.152 

0.6547 

7000 

1.68 

-0.4309 

0.900 

0.779 

0.498 

0.323 

0.128 

-0.100 

-0.274 

-0.293 

0.6766 

6900 

1.68 

-0.4086 

0.922 

0.803 

0.504 

0.317 

0.109 

-0.134 

-0.322 

-0.342 

0.6898 

6800 

1.68 

-0.3871 

0.946 

0.830 

0.510 

0.310 

0.089 

-0.168 

-0.369 

-0.391 

0.7030 

6600 

1.68 

-0.3430 

1.001 

0.891 

Z 

0.525 

=0.02 Y: 

0.296 

=0.27 M= 

0.049 

0.54 M 0 

-0.238 

-0.467 

-0.491 

0.7305 

7100 

1.49 

-0.5915 

1.393 

1.241 

0.945 

0.774 

0.598 

0.390 

0.237 

0.218 

0.5681 

7000 

1.49 

-0.5690 

1.419 

1.267 

0.950 

0.766 

0.576 

0.355 

0.188 

0.168 

0.5828 

6900 

1.49 

-0.5505 

1.444 

1.292 

0.955 

0.758 

0.557 

0.323 

0.144 

0.123 

0.5948 

7000 

1.59 

-0.4864 

1.172 

1.013 

0.696 

0.511 

0.324 

0.103 

-0.060 

-0.081 

0.6323 

6900 

1.59 

-0.4644 

1.197 

1.038 

0.701 

0.504 

0.304 

0.070 

-0.106 

0.128 

0.6441 

6800 

1.59 

-0.4436 

1.228 

1.067 

0.707 

0.497 

0.284 

0.036 

-0.153 

-0.176 

0.6580 
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Table 5. Analytical relations for the edges of the instability strip into the observational plane. 
They have the form: Color erfge =a*Mag + b *log Z + c. 


Edge 

Mag 

Color 

a 

b 

C 

FOBE 

My 

(B - V) 

-0.05±0.02 

0.05±0.()1 

0.39±0.01 



(V-/) 

-0.09±0.01 

0.03±0.()1 

0.45±0.01 


M r 

( J-K) 

-0.03±0.01 


0.17±0.01 

FRE 

My 

(B — V) 

-0.12±0.07 

0.15±0.03 

0.98±0.11 



( V-I ) 

-0.13±0.02 

0.04±0.01 

0.86±0.03 


Mjc 

(J-K) 

-0.07±0.01 


0.32±0.01 


Table 6. Coefficients of the predicted metal-dependent optical and NIR ( RIJHK ) PLZ relations 
for FU, FO and FU+FO pulsators. They take account of the entire metallicity range 
(Z=0.0001-0.02). The global (FU+FO) relations were derived by fundamentalizing the FO 

periods using the canonical relation. 


Type a 

Band b 

a c 

b A 

c e 

crj 

(Tb f 

crj 

rms g 

FU 

R 

0.24 

-1.39 

0.17 

0.02 

0.07 

0.01 

0.10 

FO 

R 

-0.34 

-2.00 

0.14 

0.07 

0.10 

0.01 

0.07 

Glob 

R 

0.25 

-1.22 

0.17 

0.02 

0.06 

0.01 

0.10 

FU 

I 

-0.07 

-1.66 

0.17 

0.02 

0.06 

0.01 

0.09 

FO 

I 

-0.65 

-2.20 

0.14 

0.06 

0.10 

0.01 

0.05 

Glob 

I 

-0.07 

-1.53 

0.17 

0.02 

0.04 

0.01 

0.09 

FU 

J 

-0.51 

-1.98 

0.17 

0.01 

0.04 

0.01 

0.06 

FO 

J 

-1.07 

-2.46 

0.15 

0.04 

0.08 

0.01 

0.04 

Glob 

J 

-0.50 

-1.90 

0.18 

0.01 

0.03 

0.01 

0.06 

FU 

H 

-0.76 

-2.24 

0.19 

0.01 

0.02 

0.01 

0.04 

FO 

H 

-1.31 

-2.70 

0.16 

0.02 

0.04 

0.01 

0.02 

Glob 

H 

-0.76 

-2.22 

0.18 

0.01 

0.02 

0.01 

0.04 

FU 

K 

-0.82 

-2.27 

0.18 

0.01 

0.02 

0.01 

0.03 

FO 

K 

-1.37 

-2.72 

0.15 

0.02 

0.04 

0.01 

0.02 

Glob 

K 

-0.82 

-2.25 

0.18 

0.01 

0.02 

0.01 

0.04 


a Pulsation mode. 
b Photometric band. 
c Zero-point of the PLZ relation. 
d Coefficient of the logarithmic period. 
Coefficient of the metallicity term. 
e Errors on coefficients. 
f Standard deviation. 
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Table 7. Coefficients of the predicted dual-band, metal-dependent, optical and NIR PWZ 
relations for FU, FO and FU+FO pulsators. They take account of the entire metallicity 
range (Z=0.00010.02). The global (FU+FO) relations were derived by fundamentalizing the 

FO periods using the canonical relation. 


Mode a 

Band b 

f 

a d 

b e 

c f 

(T ^ 
u a 

<T h S 

cr c g 

rms h 

FU 

V.B-V 

3.06 

-1.11 

-2.67 

-0.02 

0.02 

0.04 

0.01 

0.07 

FO 

V.B-V 

3.06 

-1.47 

-2.65 

0.02 

0.05 

0.09 

0.01 

0.04 

Glob 

V.B-V 

3.06 

-1.07 

-2.49 

0.01 

0.01 

0.04 

0.01 

0.08 

FU 

R,B-R 

1.68 

-1.19 

-2.74 

0.03 

0.01 

0.04 

0.01 

0.07 

FO 

R.B R 

1.68 

-1.62 

-2.82 

0.05 

0.04 

0.08 

0.01 

0.04 

Glob 

R,B-R 

1.68 

-1.17 

-2.61 

0.04 

0.01 

0.03 

0.01 

0.07 

FU 

R,V-R 

4.94 

-1.30 

-2.82 

0.09 

0.01 

0.04 

0.01 

0.06 

FO 

R.V-R 

4.94 

-1.79 

-3.02 

0.08 

0.03 

0.07 

0.01 

0.03 

Glob 

R.V-R 

4.94 

-1.29 

-2.76 

0.09 

0.01 

0.03 

0.01 

0.06 

FU 

I.B-I 

0.78 

-0.98 

-2.49 

0.11 

0.01 

0.02 

0.01 

0.04 

FO 

I.B-I 

0.78 

-1.48 

-2.77 

0.10 

0.03 

0.06 

0.01 

0.03 

Glob 

I.B-I 

0.78 

-0.97 

-2.40 

0.11 

0.01 

0.02 

0.01 

0.04 

FU 

I.VI 

1.38 

-0.94 

-2.43 

0.15 

0.01 

0.02 

0.01 

0.03 

FO 

I.VI 

1.38 

-1.49 

-2.81 

0.13 

0.02 

0.04 

0.01 

0.02 

Glob 

I.VI 

1.38 

-0.94 

-2.37 

0.15 

0.01 

0.02 

0.01 

0.04 

FU 

I.R-I 

2.31 

-0.80 

-2.28 

0.17 

0.01 

0.02 

0.01 

0.03 

FO 

I.R-I 

2.31 

-1.37 

-2.73 

0.15 

0.02 

0.04 

0.01 

0.02 

Glob 

I.R-I 

2.31 

-0.80 

-2.23 

0.17 

0.01 

0.02 

0.01 

0.04 

FU 

J.B-J 

0.27 

-0.94 

-2.36 

0.16 

0.01 

0.02 

0.01 

0.03 

FO 

J.B-J 

0.27 

-1.48 

-2.74 

0.14 

0.02 

0.05 

0.01 

0.02 

Glob 

J.B-J 

0.27 

-0.94 

-2.31 

0.16 

0.01 

0.02 

0.01 

0.04 

FU 

J.V-J 

0.40 

-0.92 

-2.33 

0.17 

0.01 

0.02 

0.01 

0.03 

FO 

J.V-J 

0.40 

-1.48 

-2.74 

0.15 

0.02 

0.04 

0.01 

0.02 

Glob 

J.V-J 

0.40 

-0.93 

-2.29 

0.17 

0.01 

0.02 

0.01 

0.03 

FU 

J.R-J 

0.52 

-0.89 

-2.29 

0.18 

0.01 

0.02 

0.01 

0.03 

FO 

J.R-J 

0.52 

-1.45 

-2.72 

0.15 

0.02 

0.05 

0.01 

0.02 

Glob 

J.R-J 

0.52 

-0.90 

-2.25 

0.18 

0.01 

0.02 

0.01 

0.04 

FU 

J.I-J 

0.96 

-0.92 

-2.29 

0.18 

0.01 

0.02 

0.01 

0.03 

FO 

J.I-J 

0.96 

-1.47 

-2.72 

0.15 

0.02 

0.05 

0.01 

0.02 

Glob 

J.I-J 

0.96 

-0.92 

-2.25 

0.18 

0.01 

0.02 

0.01 

0.04 

FU 

H.B-H 

0.16 

-1.04 

-2.50 

0.18 

0.01 

0.02 

0.01 

0.03 

FO 

H.B-H 

0.16 

-1.58 

-2.89 

0.15 

0.02 

0.03 

0.01 

0.01 

Glob 

H.B-H 

0.16 

-1.05 

-2.50 

0.17 

0.01 

0.01 

0.01 

0.03 

FU 

H.V-H 

0.22 

-1.04 

-2.49 

0.19 

0.01 

0.02 

0.01 

0.03 

FO 

H.V-H 

0.22 

-1.59 

-2.90 

0.16 

0.01 

0.03 

0.01 

0.01 

Glob 

H.V-H 

0.22 

-1.05 

-2.50 

0.18 

0.01 

0.01 

0.01 

0.03 

FU 

H.R-H 

0.27 

-1.03 

-2.47 

0.19 

0.01 

0.02 

0.01 

0.03 

FO 

H.R-H 

0.27 

-1.58 

-2.89 

0.16 

0.01 

0.03 

0.01 

0.01 

Glob 

H.R-H 

0.27 

-1.04 

-2.49 

0.19 

0.01 

0.01 

0.01 

0.03 

FU 

H.I-H 

0.44 

-1.06 

-2.50 

0.19 

0.01 

0.02 

0.01 

0.03 

FO 

H.I-H 

0.44 

-1.61 

-2.92 

0.16 

0.01 

0.03 

0.01 

0.01 

Glob 

H.I-H 

0.44 

-1.07 

-2.52 

0.19 

0.01 

0.01 

0.01 

0.03 

FU 

H.J-H 

1.68 

-1.17 

-2.67 

0.21 

0.01 

0.02 

0.01 

0.04 

FO 

H.J-H 

1.68 

-1.72 

-3.09 

0.17 

0.02 

0.04 

0.01 

0.02 

Glob 

H.J-H 

1.68 

-1.19 

-2.76 

0.20 

0.01 

0.02 

0.01 

0.04 

FU 

K,B-K 

0.10 

-1.00 

-2.44 

0.18 

0.01 

0.02 

0.01 

0.03 
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Table 7—Continued 


Mode a 

Band b 

f 

a d 


c f 

cr a g 

<T b g 

cr c g 

rms h 

FO 

K.B K 

0.10 

-1.54 

-2.84 

0.15 

0.02 

0.03 

0.01 

0.02 

Glob 

K.B-K 

0.10 

-1.01 

-2.43 

0.17 

0.01 

0.01 

0.01 

0.03 

FU 

K.V-K 

0.13 

-1.00 

-2.43 

0.18 

0.01 

0.02 

0.01 

0.03 

FO 

K.V-K 

0.13 

-1.55 

-2.84 

0.15 

0.02 

0.03 

0.01 

0.02 

Glob 

K.V-K 

0.13 

-1.01 

-2.43 

0.18 

0.01 

0.02 

0.01 

0.03 

FU 

K.R-K 

0.16 

-0.99 

-2.42 

0.19 

0.01 

0.02 

0.01 

0.03 

FO 

K.R-K 

0.16 

-1.54 

-2.84 

0.16 

0.02 

0.03 

0.01 

0.01 

Glob 

K.R-K 

0.16 

-1.00 

-2.42 

0.18 

0.01 

0.01 

0.01 

0.03 

FU 

K.I-K 

0.25 

-1.00 

-2.43 

0.19 

0.01 

0.02 

0.01 

0.03 

FO 

K.I-K 

0.25 

-1.55 

-2.85 

0.16 

0.02 

0.03 

0.01 

0.01 

Glob 

K.I-K 

0.25 

-1.01 

-2.44 

0.18 

0.01 

0.01 

0.01 

0.03 

FU 

K.J-K 

0.69 

-1.03 

-2.48 

0.19 

0.01 

0.02 

0.01 

0.03 

FO 

K.J-K 

0.69 

-1.58 

-2.89 

0.16 

0.02 

0.03 

0.01 

0.01 

Glob 

K.J-K 

0.69 

-1.05 

-2.50 

0.18 

0.01 

0.01 

0.01 

0.03 

FU 

K.H-K 

1.87 

-0.94 

-2.34 

0.18 

0.01 

0.02 

0.01 

0.03 

FO 

K.H-K 

1.87 

-1.48 

-2.75 

0.15 

0.02 

0.04 

0.01 

0.02 

Glob 

K.H-K 

1.87 

-0.94 

-2.32 

0.18 

0.01 

0.02 

0.01 

0.03 


a Pulsation mode. 
b Photometric band. 

c Coefficient of the color term. They have been estimated using the ? reddining 
law and the central wavelengths of the Landolt optical photometric system (?) and 
of the 2MASS NIR photometric system (?). The above changes also imply a total 
to selective absorption ratio Ry = 3.06. The new assumptions will be addressed in 
a forthcoming paper. 

d Zero-point of the PLZ relation. 

e Coefficient of the logarithmic period. 

Coefficient of the metallicity term. 

g Errors on coefficients. 

h Standard deviation. 
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Table 8. Same as Table [71 but for triple band PWZ relations. 


Mode a 

Band b 

e 

« d 

b c 

c f 

(T g 
u a 

o- b z 

°7 g 

rms h 

FU 

V,B-I 

1.34 

-1.01 

-2.54 

0.08 

0.01 

0.03 

0.01 

0.05 

FO 

V,B-I 

1.34 

-1.48 

-2.74 

0.08 

0.03 

0.07 

0.01 

0.03 

Glob 

V,B-I 

1.34 

-0.99 

-2.42 

0.08 

0.01 

0.02 

0.01 

0.05 

FU 

V,B R 

2.02 

-1.17 

-2.72 

0.02 

0.01 

0.04 

0.01 

0.07 

FO 

V,B-R 

2.02 

-1.58 

-2.78 

0.04 

0.04 

0.08 

0.01 

0.04 

Glob 

V,B R 

2.02 

-1.15 

-2.58 

0.03 

0.01 

0.03 

0.01 

0.07 

FU 

J.B-R 

0.57 

-0.99 

-2.44 

0.13 

0.01 

0.02 

0.01 

0.03 

FO 

J.B-R 

0.57 

-1.51 

-2.75 

0.12 

0.03 

0.05 

0.01 

0.02 

Glob 

J.B-R 

0.57 

-0.99 

-2.37 

0.13 

0.01 

0.02 

0.01 

0.04 

FU 

J.B I 

0.38 

-0.95 

-2.39 

0.14 

0.01 

0.02 

0.01 

0.03 

FO 

J.B I 

0.38 

-1.48 

-2.74 

0.13 

0.02 

0.05 

0.01 

0.02 

Glob 

J.B I 

0.38 

-0.95 

-2.33 

0.15 

0.01 

0.02 

0.01 

0.04 

FU 

J.V-I 

0.68 

-0.93 

-2.36 

0.17 

0.01 

0.02 

0.01 

0.03 

FO 

J.V-I 

0.68 

-1.48 

-2.76 

0.14 

0.02 

0.04 

0.01 

0.02 

Glob 

J.V-I 

0.68 

-0.93 

-2.31 

0.17 

0.01 

0.02 

0.01 

0.03 

FU 

H.B-R 

0.35 

-1.05 

-2.52 

0.16 

0.01 

0.02 

0.01 

0.03 

FO 

H.B R 

0.35 

-1.58 

-2.87 

0.14 

0.02 

0.04 

0.01 

0.02 

Glob 

H.B R 

0.35 

-1.06 

-2.51 

0.16 

0.01 

0.01 

0.01 

0.03 

FU 

H.B I 

0.24 

-1.03 

-2.49 

0.17 

0.01 

0.02 

0.01 

0.03 

FO 

H.B I 

0.24 

-1.57 

-2.87 

0.14 

0.02 

0.03 

0.01 

0.02 

Glob 

H.B I 

0.24 

-1.04 

-2.49 

0.16 

0.01 

0.01 

0.01 

0.03 

FU 

H.V-I 

0.42 

-1.02 

-2.48 

0.18 

0.01 

0.02 

0.01 

0.03 

FO 

H.V-I 

0.42 

-1.57 

-2.88 

0.15 

0.02 

0.03 

0.01 

0.01 

Glob 

H.V-I 

0.42 

-1.03 

-2.48 

0.18 

0.01 

0.01 

0.01 

0.03 

FU 

K.B-R 

0.23 

-1.02 

-2.46 

0.16 

0.01 

0.02 

0.01 

0.03 

FO 

K.B-R 

0.23 

-1.55 

-2.84 

0.14 

0.02 

0.04 

0.01 

0.02 

Glob 

K.B-R 

0.23 

-1.02 

-2.45 

0.16 

0.01 

0.01 

0.01 

0.03 

FU 

K.B I 

0.16 

-1.00 

-2.44 

0.17 

0.01 

0.02 

0.01 

0.03 

FO 

K.B I 

0.16 

-1.54 

-2.83 

0.15 

0.02 

0.04 

0.01 

0.02 

Glob 

K.B I 

0.16 

-1.01 

-2.43 

0.17 

0.01 

0.01 

0.01 

0.03 

FU 

K.V-I 

0.28 

-0.99 

-2.43 

0.18 

0.01 

0.02 

0.01 

0.03 

FO 

K.V-I 

0.28 

-1.54 

-2.84 

0.15 

0.02 

0.03 

0.01 

0.02 

Glob 

K.V-I 

0.28 

-1.00 

-2.42 

0.18 

0.01 

0.01 

0.01 

0.03 


a Pulsation mode. 
b Photometric bands. 
c Coefficient of the color term. 
d Zero-point of the relation. 
e Coefficient of the logarithmic period. 
f Coefficient of the metallicity term. 
g Errors on coefficients. 
h Standard deviation. 
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Table 9. Coefficients of the predicted metal-independent, dual and triple band, optical PWZ 

relations for FU, FO and FU+FO pulsators. 


Mode a 

Band b 


a d 

b e 

crj 

o~b f 

rras g 

FU 

V,B-R 

2.02 

-1.20 

-2.74 

0.01 

0.04 

0.07 

FO 

V,B-R 

2.02 

-1.74 

-3.00 

0.04 

0.10 

0.05 

Glob 

V,B-R 

2.02 

-1.20 

-2.60 

0.01 

0.03 

0.07 

FU 

V,B-V 

3.06 

-1.08 

-2.66 

0.01 

0.04 

0.07 

FO 

V.B-V 

3.06 

-1.56 

-2.76 

0.04 

0.09 

0.05 

Glob 

V,B-V 

3.06 

-1.07 

-2.49 

0.01 

0.04 

0.08 

FU 

R.B-R 

1.68 

-1.24 

-2.76 

0.01 

0.04 

0.07 

FO 

R.B-R 

1.68 

-1.79 

-3.00 

0.04 

0.10 

0.05 

Glob 

R.B-R 

1.68 

-1.24 

-2.63 

0.01 

0.03 

0.07 


a Pulsation mode. 
b Photometric bands. 
c Coefficient of the color term. 
d Zero-point of the relation. 
e Coefficient of the logarithmic period. 
f Errors on coefficients 


g Standard deviation. 



